AN  OPTIMIZATION  STUDY  OP 
BLAST  SHELTER  DEPLOYMENT 


D.ivid  L.  Mitchell 


VOLUME  III 


Appendix  H  -  BLAST  -  The  Computer  Progra 


September  1,  1966 


OCD  Contract  PS-66-1 13,  Work  Unit  1632A 


LAMBDA  REPORT  3 


AN  OPTIMIZATION  STUDY  OF  BLAST 
SHELTER  DEPLOYMENT 

by  David  L.  Mitchell 


VOLUME  III:  APPENDIX  H  -  BLAST  -  THE  COMPUTER  PROGRAM 


September  1,  1966 


PREPARED  UNDER 

INSTITUTE  FOR  DEFENSE  ANALYSES 
SUBCONTRACT  138-5 

FOR 

OFFICE  OF  CIVIL  DEFENSE 
OFFICE  OF  THE  SECRETARY  OF  THE  ARMY 
WASHINGTON,  D.C.  23310 

CONTRACT  OCD-PS-66-1 13 
WORK  UNIT  1632A 


THIS  REPORT  HAS  BEEN  REVIEWED  IN  THE  OFFICE  OF  CIVIL  DEFENSE  AND 
APPROVED  FOR  PUBLICATION.  APPROVAL  DOES  NOT  SIGNIFY  THAT  THE 
CONTENTS  NECESSARILY  REFLECT  THE  VIEWS  AND  POLICIES  OF  THE  OFFICE 
OF  CIVIL  DEFENSE. 


LAMBDA  CORPORATION 
1J01  WILSON  BOULEVARD 
ARLINGTON,  VIRGINIA  22209 


THIS  DOCUMENT  HAS  BEEN  APPROVED  FOR  PUBLIC  RELEASE  AND  SALE. 
ITS  DISTRIBUTION  IS  UNLIMITED. 

Best  Available  Copy 


TABLE  OF  CONTENTS 


Page 

I.  Introduction . .  .  .  .  1 

A.  Generation  of  Defense . 2 

B.  Evaluation  of  Defense . 6 

C.  The  Bounding  Procedure . 7 

II.  Use  of  the  Program  . . 9 

A.  Basic  Program  Operation .  9 

B.  Catalog  of  Run  Options . 12 

C.  Glossary  of  Variables  in  Common . 43 

D.  Input  Data  Sequence  Definitions . 51 

E.  Values  of  Important  Program  Constants . 62 

III.  Sample  Deck .  63 

A.  Basic  Data  Deck . 63 

B.  Shelter  Input  Deck . 72 

C.  Sample  Run  Deck  *1 . 83 

D.  Sample  Run  Deck  *2 . 92 

E.  Sample  Run  Deck  43  .  . . 98 

F.  Sample  Run  Deck  #4 . 101 

G.  Sample  Run  Deck  *5 . .  .  104 

IV.  Subroutine  Descriptions . 107 


I.  INTRODUCTION 


The  basic  tool  for  the  blast  shelter  study  is  a  computer  program,  BLAST. 

This  program  performs  three  functions  which  operate  independently,  but  are  contained 
in  the  same  program  so  that  subroutines  and  data  may  be  shared  and  so  that  any  common 
parts  will  be  consistent. 

The  first  function  of  the  program  is  to  generate  a  blast  shelter  defense  posture  of 
a  specified  cost  level  for  whatever  population  data  is  supplied,  and  then  to  calculate 
the  fatalities  which  that  posture  permits  for  a  range  of  attack  levels.  The  second  func¬ 
tion  of  the  program  is  to  compute,  for  any  cost  level  specified,  upper  and  lower  bounds 
for  the  optimum  defense  performance  at  each  attack  level  -  the  bounding  procedure. 

The  third  function  is  to  evaluate  independently  the  cell  model  by  comparing  results 
using  the  cell  model  with  results  using  a  separate  laydown  method  -  the  targeting  model. 
Since  the  targeting  model  is  used  quite  separately  ?rom  the  rest  of  the  proyam,  that  por¬ 
tion  of  the  program  is  discussed  in  Appendix  B . 

The  program  itself  has  the  same  three  divisions  with  a  driver  program  which  control 
the  selection  of  the  program  to  be  run.  If  the  driver  data  is  read  in,  parameters  are  set, 
and  then  the  portion  of  the  program  containing  the  desired  program  (called  an  overlay)  is 
read  in.  The  desired  program  is  executed  using  subroutines  in  the  main  portion  of  the 
program  plus  the  proper  overlay. 

The  next  three  sections  describe  how  to  use  the  program  in  detail.  The  first 
describes  the  options  and  parameters  to  be  set,  the  second  gives  sample  data  decks  for 
the  different  runs  and  the  data  sequences  for  all  possible  types  of  runs,  and  the  third 
describes  the  operations  of  the  subroutines. 


A.  GENERATION  OF  DEFENSE 

The  program  is  designed  to  find  a  defense  which  minimizes  fatalities, subject  to 
a  cost  constraint.  Each  possible  defense  for  each  individual  cell  is  tested  under  an 
attack  optimized  aganst  that  ceil  with  that  defense  posture;  and  it  is  this  optimum  kill 
that  is  minimized  by  the  program  for  a  given  cost.  How  the  program  does  this  is  described 
later  in  this  Appendix.  The  purpose  here  is  to  discuss  the  parameters  which  may  be  used 
to  vary  that  defense. 

Aside  from  the  population  data  and  the  shelter  and  hardness  data,  only  three 
parameters  concerning  the  defense  need  to  be  set.  The  fiist  is  TCOST,  the  cost  of  the 
nationwide  posture.  The  second  is  IFILL,  which  is  the  shelter  occupancy  or  filling 
mode.  IFILL  =  1  is  an  optimal  filling  mode;  i.e.,  fill  shelters  starting  with  the 
hardest  until  either  people  or  shelters  run  out,  and  assign  any  remaining  poonlation 
to  the  lowest  shelter  category.  IFILL  =  2  is  a  different  filling  mode,  which  assumes 
that  some  specified  fraction  of  the  population  is  alerted.  Those  alerted  fill  shelters 
in  proportion  to  the  number  of  spaces  available;  that  is,  if  the  shelters  are  underfilled, 
each  shelter  is  filled  to  the  same  percentage.* 

To  understand  the  third  parameter,  we  need  to  give  consideration  to  what  choices 
may  be  made  for  the  defense  on  each  cell . 

*A  third  possibility  for  IFILL,  for  use  with  balanced  defense,  has  no  interpretation  as 
a  filling  mode;  rather  it  was  the  most  convenient  method  for  implementing  balanced 
defense. 
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The  defense  options  are  a  number  (15  or  less  as  the  program  stands)  of  local 
shelter  mixes  available  for  use  on  each  cell .  Each  shelter  mix  defines  the  percentage 
of  the  population  in  that  cell  to  be  assigned  each  type  of  shelter.  The  entire  nationwide 
defense  consists,  then,  of  a  set  of  numbers  (from  1  to  15)  specifying  which  shelter 
mix  is  assigned  to  each  cell.  The  point  of  this  discussion  is  that  since  there  are  two 
populations  considered  in  this  model  (night  and  day),  it  must  be  specified  which  popula¬ 
tion  is  to  be  used  for  each  cell  in  determining  how  many  shelters  is  meant  by  some 
defense  option*.  The  rule  for  selection  of  the  population,  IDPOP,  to  be  used  across 
the  nation  serves  as  a  basis  for  determining  the  number  of  shelters  to  be  assigned  to 
each  cell.  Of  the  many  possible  rules,  four  are  available  in  this  program: 

IDPOP  =  1  in  each  cell,  select  nighttime  (resident)  population-called  NIGHT 
=  2  in  each  cell,  select  daytime  populatiomcalled  DAY 
=  3  in  each  cell,  select  whichever  is  larger,  night  or  dayscalled  MAX 
=  4  in  each  cell,  select  whichever  is  smaller,  night  or  dayrcalled  MIN 
Although  TCOST,  IFILL,  and  IDPOP  are  the  only  three  parameter  choices  which 
need  to  be  made  concerning  the  defense,  it  must  be  observed  that  if  Lhe  program  is  to 
minimize  fatalities  in  an  optimized  attack,  some  assumptions  need  to  be  made  about 
the  aims  of  the  optimized  attack.  Note  that  this  is  not  the  same  as  specifying  an 
entire  attack  in  advance;  this  is  merely  specifying  the  goals  to  be  achieved  in  the  attack 
optimization.  There  are  three  major  assumptions  to  be  made. 

The  first  assumption  about  the  attack  is  an  assumed  time  of  attack.  This  choice 
is  made  by  setting  the  variable  IAP0P.  The  values  for  IAPOP  are  1,2,3  and  4  with 


*For  example,  if  defense  option  11  is  70%  of  the  population  supplied  100  psi  shelters, 
and  if  cell  58  Has  100,000  people  in  the  day  and  53,000  at  night,  we  need  to  know 
whether  defense  option  11  means  assign  70,000  shelters  or  37,100  shelters. 
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the  same  definitions  as  for  IDPOP.  in  contrast  to  IDPOP,  of  course,  the  attacker 
cannot  really  choose  a  mixed  rule  for  selecting  population,  such  as  values  3  or  4;  his 
only  choice  is  to  select  a  time  of  attack,  night  or  day.  However,  these  "mixed  rules" 
for  IAPOP  generate  defenses  which  may  be  of  interest;  hence,  they  are  included  as 
possible  values.  To  create  a  defense,  one  must  set,  in  addition  to  other  parameters, 
both  IDPOP  and  IAPOP.  Selections  for  these  parameters  are  written  in  the  form 
MAX/DAY,  meaning  the  choice  for  IDPOP  is  3  and  foi  IAPOP  is  2. 

The  second  assumption  about  the  attack  concerns  the  objective  of  the  attacker . 

The  attacker  optimizes  his  payoff,  which  may  be  either  population  or  industry,  or  any 
linear  combination  of  the  two.  The  program  can  handle  up  to  10  of  these  linear  com¬ 
binations,  and  these  are  referred  to  as  the  attack  objectives,  each  of  whicn  is  defined 
by  a  population  multiplier  POPMULT,  and  an  industrial  multiplier,  VINDMULT. 

Then,  in  making  the  shelter  deployment,  the  defense  may  assume  any  one  of  those 
attack  objectives  or  he  may  assume  any  number  of  them  with  a  weight  or  confidence 
attached  to  each.  This  confidence  is  defined  in  an  array  WDAOB,  the  weight  the 
defense  assigned  to  each  attack  objective.  See  the  sample  run  decks  for  examples 
of  these  attack  objectives. 

The  third  assumption  the  defense  must  make  about  the  attack  concerns  the  attack 
level,  or  total  number  of  equivalent  megatons  expected  to  be  delivered.  A  defense 
optimized  against  a  few  megatons  is  vastly  different  from  one  optimized  against  10,000. 
There  is,  however,  a  problem  with  this  assumption.  The  basic  population  data  and 
weapon  effects  model  is  constructed  on  the  basis  of  separate  and  independent  cells. 

Only  in  this  way  can  a  defense  be  found  to  minimize  the  maximum  number  of  fatalities 
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using  the  double  Lagrange  technique;  lienee,  any  assumptions  made  about  the  attack 
level  must  be  in  a  separable  form  also.  Instead  of  assuming  an  attack  level,  the  double 
Lagrange  technique  assumes  a  lambda  to  determine  the  shelter  allocation.  This  lambda 
is  the  attacker's  marginal  or  minimum  required  fatalities  per  megaton.  The  attacker 
attacks  each  and  every  cell  in  an  amount  so  that  the  marginal  number  of  fatalities  per 
megaton  equals  lambda.  The  double  Lagrange  method  does  insure  that  the  defense 
obtained  is  best  for  whatever  attack  level  results  from  the  lambda  used. 

The  program  is  more  flexible  than  the  strict  double  Lagrange  technique.  In  order 
to  obtain  a  defense  which  is  good  against  a  range  of  attack  levels,  the  defender  is 
allowed  to  assume  a  number  of  lambdas  (DLAM's)  and  to  assign  to  each  a  weight 
(WDLAM's).  This  corresponds  tr  finding  a  defense  which  tries  to  be  good  against  a 
number  of  attack  levels,  since  each  lambda  corresponds  to  some  attack  level.  Or  the 
defense  may  choose  a  number  of  lambdas  and  try  to  minimize  the  fatalities  at  the  worst 
lambda— the  minimax  selection  of  defense. 

As  to  exactly  how  the  lambdas  are  chosen— a  small  amount  of  trial  and  error  with 
the  aid  of  an  evaluation  or  a  bounding  run  provides  sufficient  information  about  what 
lambdas  correspond  to  what  attack  levels,  and  the  exact  values  of  the  lambdas  appear 
not  to  be  important.  As  for  what  mixed  lambda  combination  (lambdas  plus  weights)  or 
mtnimax  combination  results  in  good  stabilized  defenses,  no  simple  procedure  was 
found.  It  may  be  useful  to  experiment  with  small  data  first,  and  to  apply  a  good  com¬ 
bination  to  the  large-scale  data.  Or  perhaps  a  few  trials  with  the  large-scale  data  may 
find  a  sufficiently  good  combination.  It  was  found  that  two  lambdas  for  either  the  mixed 
lambda  or  minimax  technique  seemed  to  be  sufficiently  flexible  and  much  easier  to  handle 
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than  a  larger  number  of  lambdas. 

One  last  statement  about  these  assumptions — the  justification  for  any  defense 
is  in  the  evaluation.  Any  method  of  finding  a  defense  including  writing  down  random 
numbers  is  valid  if  it  gives  desired  results  in  the  evaluation.  All  the  above  parameters 
with  the  exception  of  TC05T  and  IFILL  need  no  justification  besides  that  test,  and 
hence  one  need  not  worry,  for  instance,  if  the  attacker  really  will  have  a  mixed  attack 
objective,  if  that  assumption  gives  good  results  against  the  various  possibilities  of 
attack . 

B.  EVALUATION  OF  DEFENSE 

When  a  nationwide  defense  has  been  computed  or  entered  into  memory  from  cards, 
the  next  step  is  to  evaluate  the  defense.  For  each  of  a  number  of  attack  levels  (total 
delivered  weapons  in  terms  cf  one-megaton  equivalents),  the  program  finds  the  optimum 
attack  against  the  nationwide  defense  and  computes  the  payoff  (combination  of  population 
and  industry)  and  the  population  and  industry  killed.  Most  of  the  same  parameters 
need  to  be  set  as  for  generation  of  the  defense;  however,  the  interpretation  is  different— 
they  are  now  actual  objectives  attributed  to  the  attacker.  The  selection  of  tliese  ob¬ 
jectives  and  parameters  is  independent  of  what  was  assumed  by  the  defense. 

The  first  parameter  is  the  actual  time  of  attack,  which  is  set  in  IAPOP,  but  which 
differs  from  the  previous  use  of  IAPOP,  which  was  only  for  the  purposes  of  generating 
tlie  defense.  Only  two  realistic  choices  exist  for  actual  time  of  attack:  IAPOP  =  1  for 
nighttime  attack  and  IAPOP  -  2  for  daytime  attack.  In  order  to  evaluate  a  defense, 
both  times  of  attack  are  usually  used  in  turn. 

The  second  parameter  is  the  attacker's  objective— the  relative  weighting  of  popula¬ 
tion  and  industry  killed.  The  attack  objectives  are  defined  by  one  or  more  values  in 
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the  arrays  POPMULT  and  VINDMULT,  each  combination  corresponding  to  one  attack 
objective.  Only  the  relative  size  of  POPMULT  and  VINDMULT  matters  here,  and  not 
the  absolute  values. 

Also,  for  another  evaluation  of  a  defense  one  might  vary  IF  ILL  from  optimal  to 
fractional  filling  modes,  (See  the  previous  section.)  Other  evaluations  might  be  made 
for  softer  shelters  (via  the  PERTURB  option),  for  variations  tn  the  defmttion  of  un¬ 
sheltered  population,  or  for  different  height  of  burst  selections. 

C.  THE  BOUNDING  PROCEDURE* 

For  any  size  defense  budget,  we  would  be  interested  in  knowing  to  what  level 
the  defense  can  hold  fatalities.  When  this  minimum  is  determined  for  a  number  of  attack 
levels,  the  result  can  be  called  the  optimal  curve.  It  represents  (for  some  set  of  para¬ 
meters)  the  lowest  level  of  fatalities  the  defender  can  achieve  against  an  optimized 
attack  of  known  total  size.  We  cannot  actually  compute  the  optimal  curve,  but  we  can 
find  upper  and  lower  bounds  for  it,  and  these  can  be  made  as  close  together  as  desired 
for  some  expenditure  of  computation  time.  The  method  for  the  computation  is  described 
in  Appendix  A. 

There  are  not  many  parameters  to  set  for  the  procedure;  cost,  of  course,  is  the 
important  parameter.  IDPOP  and  IAPOP  need  to  be  set;  the  actual  time  of  attack 
will  always  be  the  same  as  the  assumed  time  of  attack.  Any  attack  objective  may  be 

*The  bounding  procedure  described  here  was  originally  conceived  as  a  method  for  checking 
solutions  generated  by  the  double  Lagrange  procedure  to  see  if  they  were  spurious.  How¬ 
ever  for  a  defense  which  is  to  be  stabilized  rather  than  optimized  for  a  single  set  of  para¬ 
meters,  a  solution  is  not  spurious  if  it  gives  desired  results.  A  standard  of  performance 
is  still  needed,  and  that  standard  is  provided  by  the  bounding  procedure. 
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selected,  but  since  the  payoff  to  defender  and  attacker  must  be  the  same,  the  only 
attack  objective  which  will  indicate  population  fatalities  is  a  pure  population  objective. 
All  other  parameters  which  need  to  be  set  have  to  do  with  program  operation  and  have 
no  larger  significance. 
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II.  USE  OF  THE  PROGRAM 


A,  BASIC  PROGRAM  OPERATION 

There  is  considerable  freedom  allowed  in  the  operation  of  the  program,  as  can 
be  seen  by  considering  Figure  1 ,  a  flow  chart  of  the  basic  structure  of  the  program. 
The  option  selector  merely  selects  the  operations  to  be  performed;  all  input,  output, 
and  computation  is  performed  in  the  option  portions  of  the  program. 

Although  all  the  options  are  shown  in  the  flow  chart  having  equal  status,  they 
vary  drastically  in  how  much  they  perform.  Some  options  merely  set  a  single  constant 
and  others  may  compute  for  an  hour  or  more.  As  each  option  is  completed,  control 
returns  to  the  option  selector  and  the  next  option  card  is  read.  Normal  termination 
is  by  reading  a  STOP  card;  however,  the  program  will  also  terminate  if  an  option 
card  is  read  which  the  option  selector  does  not  recognize. 

Because  of  the  way  the  program  is  constructed,  any  sequence  of  options  is 
allowed.  After  each  option  card,  the  data  sequence  for  that  option  is  read  and  that 
option  is  then  performed.  Hence  the  data  deck  structure  is  as  shewn  in  Figure  2: 
the  basic  data  deck,  followed  by  a  sequence  of  option  cards  each  with  its  data  deck, 
terminated  by  a  STOP  card.  The  data  sequences  for  each  option  may  be  large  (see 
for  example  the  shelter  input  deck,  which  is  the  INPUT  option)  or  the  particular  option 
may  not  read  any  cards  at  all.  The  data  sequences  for  all  the  options  are  specified 
in  Section  II  D  of  this  Appendix. 

The  typical  mode  of  operation  is  as  shown  in  Figure  3.  This  mode  of 
operation  is  simply  a  special  case  of  the  pattern  of  Figure  2.  The  shelter  input 
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Figure  1 .  Flow  Chart  of  Program  BLAST 
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deck  is  an  INPUT  card  followed  by  an  INPUT  option  data  sequence  (an  example  of 
a  shelter  input  deck  is  given  in  the  Sample  Data  Deck  section).  The  run  deck  is  a 
sequence  of  option  cards,  each  followed  by  its  own  option  data  deck  (examples  of  run 
decks  are  shown  in  the  Sample  Data  Deck  section).  And,  as  always,  the  STOP  card 
terminates  the  deck. 
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As  can  be  seen,  it  is  the  options  and  the  option  cards  which  control  whether 
the  program  performs  large  computations  or  small  (or  even  no)  computations.  Hence 
the  next  two  sections  define  what  the  options  do  and  what  data  is  required.  The  point 
of  view  will  be  strictly  as  a  user;  the  subroutines  are  discussed  in  Section  IV. 

B.  CATALOG  OF  RUN  OPTIONS 
1 .  Primary  Computation  Options 
Run  Type  SIMPLE 

This  is  tie  basic  option  of  the  program  —  it  allows  one  to  generate  a 
defense  and  to  evaluate  it.  There  are  many  parameters  to  be  set  before  entering 
SIMPLE  where  the  defense  is  generated.  Then  the  defense  will  be  generated  in 
SIMPLE  on  the  basis  of  whatever  values  were  chosen  for  the  parameters.  If  the 
program  is  directed  to  go  ahead  and  evaluate  the  defense,  the  same  parameters  wili 
be  used  for  the  evaluation.  If  the  defense  is  not  evaluated  immediately,  after 
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generating  the  defense  in  the  SIMPLE  option  control  returns  to  the  option  selector, 
where  the  parameters  may  be  reset  and  an  evaluation  performed. 

Parameters  which  are  to  be  set  before  eriering  SIMPLE  are  the  cell 
population  data,  the  population  bases  IDPOP  and  lAPOP,  the  shelter  data,  the 
shelter  filling  mode  ((FILL),  and  the  lethal  areas,  heights  of  burst,  the  attack  weapon 
density  quanta  ,  and  the  attack  leve's  for  the  evaluation  portion  of  SIMPLE.  The 
population  data  is  set  in  the  basic  data  deck  or  in  run  type  LIVEDATA;  run  type 
OPTIONS  is  used  to  set  IDPOP,  lAPOP,  and  IFILL;  the  shelter  data,  lethal  areas, 
and  heights  of  burst  are  set  in  the  basic  data  deck  or  in  run  type  INPUT;  and  the 
attack  density  quanta  and  attack  levels  are  set  in  the  basic  data  deck.  In  addition, 
two  decisions  concerning  the  generation  of  the  defense  need  to  be  made  before  enter¬ 
ing  SIMPLE — first,  if  a  punch  output  of  the  defense  is  desired  option  PUNCH  should 
be  entered,  and  second,  if  minmax  lambda  stabilization  is  desired  option  MINIMAX 
should  be  enter.-'" 


In  run  type  SIMPLE,  the  following  parameters  are  set  explicitly; 


TCOST 

NDLAM 

DLAM's 

NAOB 

POPMULT's 

VINDMULT's 


Shelter  program  cost 

Number  of  defense  lambdas 

Lambdas  assumed  by  the  defense 

Number  of  attack  objectives 

Population  multipliers  of  the  attack  objectives 

industry  multipliers  of  the  attack  objectives 


Recall  that  in  the  defense  optimization  every  defense  is  attacked  optimally;  hence 
these  attack  parameters  must  be  set. 
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For  description  of  precisely  how  the  above  parameters  may  be  input  see 
the  Data  Sequence  Definition  section.  Since  no  weights  on  the  attack  objectives 
(WDAOB's)  have  been  read,  the  program  uses  whatever  has  been  stored  previously. 
(As  set  initially  in  the  program  constants,  WOAOB  places  the  entire  weight  for  the 
defense  on  the  first  attack  objective;  i.e.,  the  defense  assumes  the  first  attack 
objective  unless  otherwise  set  —  subsequently  in  SIMPLE  or  in  run  type  ATTOBJ.) 

At  this  point  the  program  will  generate  the  defense  and,  if  it  has  been  called 
for  by  option  PSJMCH,  it  will  retain  the  defense  on  a  punched  card  deck. 

The  defense  is  now  fixed  and  fully  defined  and  may,  for  instance,  be  plotted  on 
maps.  The  defense  is  defined  by  one  of  the  15  defense  options  for  every  cell  plus 
the  split  cell  information  for  the  one  cell  which  must  be  split  to  achieve  exact  cost 
closing,  together  with  iDPOP  which  defines  the  population  basis  to  determine  the 
number  of  shelters  in  each  cell.  The  chosen  defenses  are  stored  in  the  memory  in 
packed  form  along  with  the  cell  population  data.  The  next  step  is  to  evaluate  the 
defense. 

The  next  card  to  be  read  by  option  SIMPLE  determines  whether  the  defense 
will  be  evaluated  immediately  or  control  will  be  returned  to  the  option  selector. 

If  control  is  returned  to  the  option  selector  then  the  defense  which  is  defined  in 
memory  may  be  evaluated  later  using  oDtion  EVAL,  or  it  may  be  analyzed  using 
OUTPUT,  or  used  for  time  phasing  studies  in  TiMEFAZE,  or  the  run  may  simply 
be  terminated.  If  the  defense  is  to  be  evaluated  immediately,  the  procp*am  will 
evaluate  the  defense  for  each  of  the  NAOB  attack  objectives,  printing  the  result 
after  each  evaluation. 
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The  evaluation  consists  of  optimizing  the  offense  using  a  one-sided 
Lagrange  multiplier  maximization/  adjusting  the  lambdas  in  the  maximization  until 
sufficiently  accurate  calculation  is  made  of  the  fatalities  for  each  of  several  attack 
levels.  The  offense  parameters  required  are  the  density  quanta,  the  attack  levels,  and 
the  attack  objectives,  and  these  have  already  been  set. 

After  each  evaluation  a  print  like  Figure  4  is  made,  identifying  the 
run  recording  the  parameters  and  giving  the  results  of  the  evaluation.  The  first 
line  of  Figure  4  identifies  this  as  a  run  type  SIMPLE.  The  key  number  in  the 
next  line  is  assigned  during  the  time  the  punched  card  deck  is  being  made.  The 
same  number  is  assigned  to  both  the  printout  and  the  punch  deck  so  that  there  is 
no  confusion  about  matching  the  two*. 

The  next  block  defines  the  parameters  used  to  determine  the  defense — the 
population  bases  IDPOP  and  IAPOP,  the  shelter  filling  mode,  the  cost  for  the  shelter 
program  (TCOST),  and  the  mu  found  which  achieves  that  cost.  Next  the  lambdas 

assumed  by  the  defense  are  given  together  with  either  the  weights  for  the  mixed 
lambda  case  or  the  statement  that  mirimax  stabilization  was  used.  The  last  parameter 
of  the  defense  is  a  listing  of  all  attack  objectives  which  have  nonzero  weights 
(WDAOB),  i.e.,  all  attack  objectives  considered  by  the  defense. 

The  next  block  defines  the  parameters  used  by  the  offense— the  attack 
density  quanta  and  the  attack  objective  used  by  the  offense  for  this  evaluation. 

Next  the  payoffs  or  fatalities  are  given  for  each  attack  level.  The  attack  level 
is  given  in  equivalent  megatons.  The  payoff  killed  is  the  fraction  of  the  total 

*  The  number  assigned  is  the  first  five  digits  of  the  fraction  of  the  split  cell 
which  is  assigned  to  the  primary  defense  for  that  cell.  This  number  is  unique  to 
each  run. 
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payoff  killed  by  the  attacker,  or  the  sum  of  the  population  killed  and  the  industry 
killed,  weighted  by  the  population  and  industry  multipliers  used  by  the  offense. 
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Figure  4.  Sample  of  SIMPLE  Output 
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In  this  case  since  the  attack  objective  is  pure  population,  the  payoff  killed  is 
the  fraction  of  total  population  killed.  The  population  killed  is  in  millions  for  a 
total  of  some  209,000 ,000;  the  industry  killed  is  in  billions  of  dollars  of 
manufacturing  value  added  for  a  total  of  $97,000,000,000. 

The  lambda  column  gives  the  lambdas  found  for  each  attack  level  in 
the  evaluation  with  no  scale  factor.  This  gives  the  marginal  return  in  payoff 
for  an  additional  megaton  of  attack.  Thus,  at  attack  level  of  5000  megatons  an 
additional  megaton  adds  approximately  7500  fatalities. 

The  final  block  of  the  SIMPLE  print  indicates  what  happens  in  the  tail 
cells— the  low  density  aggregated  cells.  For  each  cell  the  defense  option  is 
given  together  with  the  attack  quantum  for  each  attack  level  considered  in  the 
evaluation.  For  example,  in  the  figure,  tail  cell  number  3128  was  undefended 
(defense  option  1)  and  for  attack  level  5000  it  received  .05  megatons  per  square 
mile  (attack  quantum  3). 

After  making  this  print,  the  evaluation  is  repeated  for  the  remaining 
attack  objectives,  and  a  similar  printout  is  made  for  each.  After  that  a  card  is 
read  to  see  if  new  WDAOB's  are  desired.  If  so,  then  the  WDAOB's  are  read,  and  a  new 
defense  is  generated  and  evaluated  for  all  attack  objectives.  Because  of  running 
time,  this  last  option  was  not  used  on  the  large-scale  data  runs  in  this  study.  When 
option  SIMPLE  terminates,  control  returns  to  the  option  selector. 

There  are  no  error  exits  or  messages  for  option  SIMPLE;  only  two 
abnormalities  may  occur:  if  the  defense  is  unable  to  spend  all  its  target  cost  for 
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the  initial  ma's,  then  it  saturates  at  some  cost  and  evaluates  the  defense 
for  that  cost.  The  achieved  cost  is  printed  along  with  a  statement  that  the 
defense  saturated  in  the  block  of  the  SIMPLE  print  which  defines  the  defense 
parameters.  Secondly,  if  in  the  evaluation  the  offense  is  not  able  to  achieve 
the  highest  attack  levels  with  lambdas  greater  than  1.0,  then  the  number  of 
attack  levels  is  reduced  and  the  highest  attack  level  is  dropped.  No  message 
occurs,  but  this  truncation  is  obvious  in  the  printout. 

When  the  SIMPLE  option  card  is  read,  control  is  transferred  to  a 
portion  of  a  driver  program  called  GVLYDAVE,  the  first  overlay  of  the  program. 

At  that  point  subroutines  in  the  main  portion  of  the  program  and  in  the  first  overlay 
are  available  for  use.  After  the  data  is  read  in  GVLYDAVE,  the  program  then  calls 
DEFOPTG  to  generate  the  defense,  DEFPUNCH  if  the  punch  output  is  desired, 

EVAL  to  perform  the  evaluations,  and  ANSPRIN  to  print  the  results  of  each 
evaluation.  For  discussion  of  the  methodology  of  these  subroutines  refer  to  the 
section  on  each  subroutine. 

Run  Type  EVAL 

This  is  uie  option  which  evaluates  any  defense  stored  in  the  cell  data  in  memory,  no 
matter  whether  computed  or  input.  This  enables  one  to  evaluate  a  defense  for  new  attack 
objectives,  time  of  attack,  or  filling  mode. 

EVAL  itself  reads  no  data  and  sets  no  parameters.  The  entire  defense,  all 
objectives,  and  all  other  parameters  must  be  set  before  entering  option  EVAL.  (Because 
of  program  construction  these  parameters  will  always  be  set  to  some  value,  but  they 
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must  be  set  to  desired  values).  The  following  parameters  must  be  properly  set  before 
EVAL: 

1.  The  cell  population  data  together  with  the  defense  to  be  evaluated  (JDEF 
defined  for  every  cell) 

2 .  The  shelter  data ,  lethal  areas  ,  heights  of  burst 

3.  iDPOP  is  presumably  already  set  to  the  proper  basis 

4.  Shelter  filling  mode  <1  FILL) 

5.  Attack  weapon  density  quanta 

6.  Attack  levels 

7.  Time  of  attack  (IAPOP) 

8.  Attack  objective  (POPMULT,  VINDMULT  ,  and  IA0B) 

With  the  above  parameters  set,  when  the  program  reads  the  EVAL  option  card, 
optimized  defenses  will  be  generated  for  each  attack  level,  and  the  result  will  be  printed 
out.  As  in  SIMPLE,  the  procedure  is  the  Lagrange  multiplier  method.  When  the  defense 
is  evaluated,  the  result  is  printed  in  the  same  form  as  the  SIMPLE  output  (Figure  4), 
except  that  the  key  number  :s  not  assigned.  The  parameters  for  the  defense  are  whatever 
remain  for  tliose  variables,  and  these  may  o»  may  not  be  meaningful .  The  time  of  attack 
(IAPOP)  is  indicated  in  the  second  line  of  the  first  block.  The  offense  parameters  are 
those  which  have  been  used  for  the  optimization,  and  the  results  are  printed  as  before. 
Option  EVAL  terminates  after  making  the  printout,  and  control  returns  to  the  option 
selector. 

There  are  no  error  exits  or  messages  from  option  EVAL;  as  in  SIMPLE,  if  the 
attack  cannot  achieve  the  highest  attack  level,  that  attack  level  is  dropped.  The 
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EVAL  option  is  performed  from  OVLYDAVE  and  uses  only  two  subroutines,  EVAL 
and  ANSPRIN. 

Run  Type  BALDEF 

This  option  enables  one  to  compute  and  evaluate  "balanced11  defenses  for 
given  cost  levels.  The  procedure  uses  both  an  upper  hardness  limitation  and  a  lower 

it 

hardness  cutoff  . 

Before  entering  BALDEF  one  should  set  the  population  data,  IDPOP,  and 
1APOP;  shelter  data  may  be  entered  in  option  BAL  >EF  or  previously.  The  data 
sequence  section  indicates  the  parameters  to  be  set  to  define  the  balanced  defense 
shelter  options;  the  variables  have  the  same  meaning  as  in  the  INPUT  sequence. 
After  the  initial  shelter  data  are  set,  one  reads  in  the  costs  and  attack  levels  to  be 
evaluated.  The  final  parameters  are  PSIMIN,  the  lower  hardness  cutoff,  PSIMAX, 
the  upper  hardness  aitoff,  and  WMIN,  the  lethal  area  corresponding  to  WM!N. 

All  shelter  cost  data  is  computed  in  subroutine  BDCELL  using  an  analytic 
approximation  to  the  cost  data,  and  hence  any  changes  must  be  made  to  (he  function. 
The  function  used  for  FINAL  shelters  is: 

Cost  per  person  =  112  and  13.5 VpSI  ; 
the  degree  of  fit  is  shown  in  the  section  on  balanced  defense. 


Lambda  Research  Memorandum  2  indicates  the  effect  of  cutoffs  and  has  a  general 
discussion  of  balanced  defense. 
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The  output  after  the  evaluation  has  the  same  format  as  the  SIMPLE  printout,  and 
all  the  quantities  have  the  same  meaning.  The  program  then  prints  the  computed  lambda 
value;  this  lambda  gives  the  hardness  of  any  cell  according  to  the  equation 

Lethal  area  =  _ h _  . 

Population  Density 

After  the  balanced  defense  is  computed,  it  can  be  evaluated  changing  any  parameters 
desired.  The  exception  is  that  the  filling  mode  must  remain  at  the  artificial  value  of 
IFILL  =  4, 

This  option  is  performed  in  OVLYDAVE  in  subroutine  BALDEF2,  and  the 
evaluation  is  performed  in  EVAL. 

Run  Type  BOUND 

This  option,  which  computes  upper  and  lower  bounds  for  the  optimal  curve  as 
a  function  of  attack  levels,  is  the  implementation  of  the  boundary  procedure  described 
in  the  methodology  section  of  the  main  report.  For  any  specified  cost  the  bounds  can 
be  generated  as  close  to  each  other  as  desired;  the  payment  for  the  accuracy  is  in  com¬ 
putation  time. 

Before  entering  BOUND,  the  two  population  bases  IDPOP  and  IAP0P  should  be 
set.  Here  the  actual  attack  time  is  the  same  as  the  assumed  attack  time.  The  attack 
objective  may  be  set,  but  pure  population  is  the  only  attack  objective  which  gives  the 
answer  in  population  fatalities;  the  attack  objective  will  already  normally  be  set  to  pure 
population. 

There  are  several  modes  of  operation  available  with  run  type  BOUND;  however, 
only  three  were  found  at  all  useful;  the  others  are  described  briefly  in  the  data  sequence 
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definition  section. 

All  the  modes  require  the  following  data: 

MODESW 

This  is  the  switch  which  selects  the  mode  of  operation. 

FMUs 

20  mus  are  usee!  for  the  boundary  procedure.  These  mus  should 
be  chosen  carefully  so  that  resultant  costs  are  close  to  the  desired 
budget  level . 

TCOST 

This  is  the  target  cost  foe  deployments. 

HSAT 

This  is  the  total  possible  payoff  ffor  pure  population  objective,  it 
is  total  population) . 

One  of  the  keys  to  getting  good  accuracy  from  this  procedure  is  to  choose  the  mus 
properly.  If  most  of  the  20  mus  can  be  selected  in  a  region  where  they  will  give  costs 
close  to  TC03T,  then  the  cost  closing  will  be  quite  good  and  the  bounds  will  be  good. 

As  many  as  possible  should  be  concentrated  in  the  region  of  interest;  however,  some  mus 
should  be  available  at  extreme  values  so  the  cost  will  always  be  bracketed.  Previous 
bound  runs  and  defenses  generated  provide  information  for  the  mu  selection. 

With  the  above  parameters  required  for  all  modes,  we  now  discuss  them  individually. 

Mode  1 .  This  mode  finds  upper  and  lower  bounds  for  the  optimal  curve  by  sweeping 
through  all  lambdas,  starting  from  some  point  in  the  middle  for  efficiency.  The  procedure 
continues  until  it  has  swept  through  all  lambdas  From  10.0  to  1,000.000.0.  Two 
additional  parameters  are  required: 

DLAM(l)  -  This  is  the  initial  lambda  used  to  compute  the  initial  point  and 
begin  the  sweep. 

EPS  -  This  is  a  parameter  controlling  the  separation  between  the  lambdas 
chosen,  and  hence  controls  computation  time  and  accuracy. 
Approximately,  it  is  the  distance  between  upper  and  lower  bounds. 

After  reading  the  data,  the  program  computes  the  initial  point  from  the  initial  lambda. 
It  then  computes  20  lambdas  and  does  one  pass  through  ah  the  cells,  finding  solutions  of 
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the  Lagrangian  tliose  lambdas  with  defenses  of  cost  approximately  TCOST.  The  lower  and 
upper  bound  theorems  are  applied  to  the  solutions,  and  the  bounds  are  computed.  A  new 
set  of  lambdas  is  chosen  to  extend  the  range,  and  the  process  continues  until  the  bounds 
are  generated. 

Figure  5  shows  the  initial  point  used  for  a  bound  run,  the  payoff,  weapon  expendi¬ 
ture,  cost,  lambda,  and  mu.  Figure  6  shows  the  output  after  the  second  pass  through  the 
cells.  The  20  lambdas  are  shown,  together  with  the  mus  giving  the  closest  cost,  and  the 
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Figure  6.  Lagrangian  Solutions  for  Pass  Number  2 
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payoff,  weapon  expenditure,  cost,  TCOST,  and  for  each  bounding  line  the  intercept 
and  saturation  point  (H°  and  A°  in  the  methodology  section).  After  each  pass  the 
bounding  lines  are  computed  and  the  lower  bound  lines  are  punched  onto  cards.  For 
example,  for  lambdas  number  1  and  2  in  Figure  6,  the  following  card  was  punched: 


The  first  number  is  the  intercept  of  the  lower  bound  line  for  zero  attack  level  (in  millions); 
the  second  is  the  slope  of  the  line  (in  thousands).  The  third  number  is  TCOST,  the 
fourth  and  fifth  are  the  two  lambdas  used  (in  thousands).  Now,  this  bounding  line  is  the 
lower  bound  for  points  on  the  optimal  curve  with  slope  between  the  two  lambdas;  the  entire 
lower  bound  is  the  lower  envelope  of  all  such  lines.  We  will  use  this  fact  when  it  comes 
to  fill  in  lambdas  (mode  5).  And  lastly,  the  upper  and  lower  bounds  generated  (values 
In  millions)  are  shown  in  Figure  7  . 

The  only  error  return  occurs  when  the  defense  saturates  in  the  computation  of  the 
initial  point.  An  error  message  is  printed  and  the  run  terminates. 
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Figure  7.  Upper  and  Lower  Bounds  for  Lagrangian  Solutions 

Mode  2  is  exactly  the  same  as  Mode  1,  except  that  the  initial  point  is  read  in,  not 
computed.  This  point  is  used  only  for  lambda  selection  and  not  in  the  bound  computation 
itself,  so  if  a  Lagrangian  solution  is  approximately  known,  it  would  be  useful  to  use  this 
mode.  The  following  additional  parameters  are  required: 


BKILLC2) 

- 

Initial  payoff 

BWEP(2) 

- 

Initial  weapon  expenditure 

BCOST(2) 

- 

Initial  cost  (this  replaces  TCOST) 

BLAM(2) 

- 

Initial  lambda  (this  replaces  DLAM(l)) 

BMU(2) 

_ 

Initial  mu 

Outputs  are  exactly  as  in  Mode  1 . 

Mode  5.  It  may  often  be  useful  to  specify  lambdas  for  the  bound  generation,  rather 
than  having  them  computed.  This  mode  performs  one  pass  with  lambdas  which  are  input 
and  computes  the  bound  from  solutions  from  those  lambdas. 

The  chief  use  for  this  mode  is  to  improve  the  accuracy  of  a  Mode  1  or  Mode  2  run. 
To  do  this,  suppose  that  for  some  portion  of  a  run,  the  error  is  unacceptably  high.  The 
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procedure  is  as  follows: 


1.  Select  two  lambdas  from  tlie  previous  run  which  define  the  inaccurate  section. 

2.  Select  3  8  lambdas  spread  between  those  two  limits.  The  20  lambdas  thus 
found  will  be  the  lambdas  (DLAM's)  for  the  new  run. 

3.  From  the  printout  of  the  mus  chosen  for  lambdas  in  the  bracket,  select  20  new 
mus.  ft  should  be  possible  to  have  very  good  cost  matching  for  these  closely 
chosen  mus . 

4.  Perform  mode  5  run 

5.  Upper  bound  at  each  point  is  now  the  lower  of  the  upper  bounds  from  the  two  mus. 

6.  To  find  the  lower  bound,  recall  that  we  must  have  bounding  lines  which  represent 
the  entire  lambda  range  ,  Thus,  we  discard  those  output  cards  from  the  first 

run  which  cover  the  lambda  range  between  the  two  limit  lambdas  chosen  above, 
and  we  substitute  the  output  cards  from  the  second  run  (which  cover  the  same 
range).  The  composite  deck  still  covers  the  entire  range,  but  the  accuracy 
will  be  improved  .  The  lower  bound  is  then  the  lower  envelope  of  the  bounding 
lines. 

The  bounding  procedure  is  performed  in  OVLAY3  with  subroutine  PPEBND  as  a 
monitor  and  input  program.  The  lambda  selection  is  done  in  PICKLAM,  the  Lagrangian 
solutions  are  found  in  DEFOPTG,  and  the  bounds  are  computed  in  BNDSET. 

Run  Type  STATCOST 

This  option  implements  the  generation  of  defenses  under  the  constraint  of 
allocating  the  defense  budget  to  states  by  population  and  optimizing  each  state  individually 
for  its  portion  of  the  budget.  Ail  the  same  parameters  are  set  in  option  STATCOST  as 
in  SIMPLE;  in  fact,  the  two  data  sequences  are  identical  (they  are  read  from  the  same  part 
of  the  program).  The  only  difference  is  that  instead  of  a  single  optimization  in  DEFOPTG 
in  the  SIMPLE  option,  the  STATCOST  option  does  the  optimization  in  subroutine 

*  It  is  possible  to  use  all  Lagrangian  solutions  in  the  critical  range  to  find  new  bounding 
lines  using  the  theorem,  but  the  above  procedure  is  simpler  and  quite  adequate. 
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STATCOST;  however,  subroutine  STATCOST  also  does  evaluations,  and  that  necessitates 
changes  to  the  run  deck  if  extra  evaluations  are  not  to  be  performed.  First,  however,  we 
will  see  what  STATCOST  does. 

The  data  input  for  STATCOST  is  precisely  the  same  as  option  SIMPLE  —  all  the 
variables  have  the  same  interpretation  with  the  exception  that  TCOST  is  now  to  be 
divided  among  the  states  and  each  state  optimized  individually.  Figure  8  displays  the 
first  portion  of  ihe  STATCOST  print,  and  it  shows,  for  instance,  that  Connecticut's 
share  of  the  totai  budget  is  $160  million.  Connecticut,  and  every  other  state,  is  then 
optimized  in  DEFOPTG  for  the  allotted  amount  of  money.  For  each  state  the  optimization 
is  performed  using  those  cells  which  have  been  explicitly  stored  in  memory  (and  hence 
meet  the  population  density  requirement). 

Now  for  the  mus  used  in  the  program,  it  may  not  be  possible  to  spend  the  entire 
amount  allocated  to  a  state  on  the  cells  which  represent  that  state*.  This  is  indicated  in 
Maine,  for  instance.  When  that  happens,  a  message  is  printed,  and  the  defense  is  set  for 

whatever  amount  can  be  spent.  This  means  that  there  will  be  some  money  unspent  out  of 
the  initial  budgeted  amount.  What  can  be,  and  has  been  done,  is  to  redo  the  calculation 
using  a  slightly  higher  initial  cost  to  bring  the  total  amount  spent  to  the  desired  amount. 
This  added  amount  to  be  budgeted  can  be  found  as  the  unspent  money  multiplied  by  the 
ratio  of  the  initial  total  budget  to  the  money  allocated  to  the  unsaturated  states.  This 
will  then  give  the  proper  amount  spent  on  the  second  pass  if  no  new  states  saturate.  As 
shown  in  the  figure,  it  took  a  budget  of  10.38  billion  to  achieve  10  billion  spent. 

"fc 

This  saturation  is  with  respect  to  the  mus  and  in  consideration  of  the  attack  which  is 
assumed;  i  .e. ,  if  a  state  is  attacked  very  lightly  even  without  shelters,  it  may  not  be 
possible  for  it  to  spend  its  budgeted  amount. 
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Figure  8.  Sample  of  STATCOST  Output 


After  finding  the  defenses  for  the  51  states,  subroutine  STATCOST  also  punches 
out  the  defense  (with  no  tail  cells)  and  then  evaluates  the  defenses  two  ways,  both  with 
IAOB  =  1.  The  first  way  is  with  no  tail  cells;  i.e. ,  an  evaluation  of  the  cells  stored 
in  memory.  The  second  way  is  with  the  tail  cells  included,  but  undefended.  This 
second  way  is  what  is  useful  for  comparison  with  unconstrained  defenses.  Answers  are 
printed  for  both  evaluations  in  exactly  the  same  form  as  the  SIMPLE  output.  After  this, 
the  evaluations  normally  performed  ir  SIMPLE  will  be  performed  if  the  NOEVAL  card  is 
not  used.  Also,  an  additional  punch  deck  may  be  obtained  with  the  PUNCH  option. 
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Because  of  the  way  the  evaluations  are  performed,  certain  modifications  must  be 
made  to  the  run  deck  if  only  pure  population  evaluations  are  desired.  A  sample  is  as 
follows: 


LIVEDATA 

JO. 375653 

OPTIONS 

2 

IAPOP 

2 

78850.0  7885.0 

ATTOBJ 

2 

2 

0.7  0.3 

1.0  0.5 

2 

0.0  1.077 

1.0  0.5 

0.0  1.0 

0.0  1.077 

PUNCH 

NOEVAL 

STATCOST 

STOP 

This  sequence  will  generate  the  defense,  punch  two  output  decks  (one  with 
tail  cells,  one  without),  and  do  two  evaluations  (one  with  tail  cells,  one  without). 

There  are  no  error  exits  or  messages  in  STATCOST.  STATCOST  is  executed 
in  OVLYDAVE,  and  after  reading  data,  subroutine  STATCOST  is  called,  which  then 
makes  use  of  DEFOPTG,  DEFPUNCH,  EVAL,  and  ANSPRIN. 

Run  Type  TIMEFAZE 

This  option  enables  one  to  find  a  defense  which  is  time-phased  with  a  larger 
defense.  Before  entering  TIMEFAZE  the  larger  defense  should  be  entered  into  memory 
by  computation  or  by  DEFSPEC .  Calling  option  TIMEFAZE  then  sets  that  old  defense  up 
for  time-phasing.  Next  the  routine  determines  which  combinations  of  defenses  are  excluded, 
as  is  shown  in  Figure  9  (the  l's  are  excluded  combinations).  To  find  the  new  defense  one 
then  uses  option  SIMPLE  in  the  normal  manner;  those  defenses  excluded  in  the  table  will 
not  be  allowed  when  the  new  defense  is  optimized.  This  routine  requires  no  data;  it  is 
executed  in  subroutine  TIMEFAZE,  which  prepares  DCONSTR  for  proper  functioning. 
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Figure  9.  Output  for  Run  Type  TIMEFAZE 

Run  Type  UNSHEL 


This  is  the  option  which  enables  one  to  compute  unsheltered  population  fatalities 
(or  in  fallout  shelters  only,  or  whatever  is  used  for  defense  option  number  1).  This 
option  is  merely  a  special  case  of  the  EVAL  option,  and  almost  all  parameters  which 
are  set  there  for  the  attacker  should  be  set  for  UNSHEL  (the  exception  is  IFILL,  which, 
of  course,  has  no  significance  here.)  No  data  is  required  in  UNSHEL;  as  in  EVAL,  all 
parameters  must  be  set  previously. 

The  procedure  used  for  UNSHEL  is  to  set  all  cells  to  JDEF  =  1  (normally 
unsheltered  defense  option)  and  to  evaluate  in  subroutine  EVAL.  The  print  (from 
ANSPRIN)  is  exactly  as  the  SIMPLE  printout,  except  that  much  of  the  information  does 
not  apply.  Then  a  correlation  printout  is  given  which  is  like  that  in  option  OUTPUT, 
and  option  UNSHEL  terminates. 
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2.  Input  «,nd  Parameter  Options 

Run  Type  ALERT 

This  option  enables  use  of  the  alternate  filling  mode  IPILL  =  2.  This  mode 
enables  one  to  specify  a  fraction  of  the  population  alerted.  Only  one  parameter  is 
read--PALERT  ,  the  fraction  alerted.  Then  the  IFILL  parameter  is  set  to  2. 

Run  Type  ATTOBJ 

This  option  enables  one  to  change  the  attack  objective.  Four  variables  are 
set-'NAOB,  the  number  of  attack  objectives;  POPMULT's,  the  population  multipliers 
of  the  different  objectives;  VINDMULT's,  the  industry  multipliers;  and  WDAOB's,  the 
weight  the  defense  assigns  to  the  attack  objectives. 

Run  Type  DEFSPEC 

This  option,  for  all  its  input  data  sequence  complexity,  is  really  just  an  input 
procedure  enabling  one  to  specify  a  defense  for  evaluation  or  analysis  .  An  input  deck 
may  be  made  up  according  to  the  definition  in  the  next  section  for  any  defense  desired 
whatever;  however,  the  primary  method  of  use  is  for  reading  in  the  defenses  saved  on 
cards  via  the  PUNCH  option. 

The  input  sequence  for  DEFSPEC  includes  almost  all  parameters  needed  to  define 
a  defense  and  evaluate  it.  One  that  is  not  defined  is  IA0B,  the  index  of  the  attack  objec¬ 
tive.  ft  may  be  set  to  a  desired  value  before  using  DEFSPEC;  if  not  set,  it  normally  has 
a  value  of  1. 

The  last  card  of  the  input  sequence  gives  one  a  choice  of  either  evaluating  the 
defense  immediately  using  ail  parameters  as  set  or  not  evaluating.  in  the  latter  case  one 
may  evaluate  after  changing  parameters ,  or  one  may  have  other  purposes.  In  either  case, 
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the  printout  has  the  same  form  as  the  ones  in  INPUT  and  SIMPLE ,  and  the  purpose  is  to 
present  the  parameters  for  the  defense . 

Previous  to  using  DEFSPEC  the  ceil  data  will  presumably  have  been  entered.  If 
the  number  of  celis  in  the  data  is  not  equal  to  the  number  defined  in  the  defense,  an  error 
message  is  printed,  but  the  program  continues.  If  the  number  of  defense  -s  too  few,  all 
remaining  cells  are  set  to  unsheltered  posture. 

DEFSPEC  is  expected  in  OVLYDAVE  using  subroutine  DEFSPEC  and,  if  evalua¬ 
tion  is  desired,  subroutine  EVAL. 

Run  Type  INPUT 

This  option  is  the  primary  means  of  entering  shelter  data  into  the  program.  It 
defines  lethal  areas,  shelter  data,  and  defense  options,  as  described  in  the  data  sequence 
section  and  the  sample  deck  section. 

The  printout  which  indicates  the  data  that  has  been  entered  is  in  Figure  10  .  This 
output  corresponds  to  the  data  in  the  sample  shelter  input  deck  for  the  FINAL  shelters. 
FINAL-2  is  simply  a  name  for  the  defense  data,  as  are  NFSS ,  10  PSI ,  30  PSI ,  etc . 
names  for  the  shelters  defined  in  the  data. 

INPUT  is  executed  from  OVLYDAVE  and  is  performed  in  subroutine  DEFPRNT, 

which  uses  SHELINP. 

Run  Type  LIVEDATA 

This  is  the  option  which  enters  the  population-industrial  data  for 
the  2  sq.nm.  cells  into  memory.  The  data  is  read  from  a  data  tape  whose 
format  will  be  described.  No  parameters  need  be  set  before  using 
LIVEDATA;  they  may  be  set  before  or  after.  However,  option  LIVEDATA 
destroys  any  defense  stored  in  memory. 
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Figure  10.  Printout  of  Data  Entered  through  Run  Type  INPUT 
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The  data  file  must  be  a  single-reel,  binary  buffered  file  with  900 
words  per  physical  record.  The  file  is  unlabeled  and  the  terminal  record 
is  to  be  filled  with  the  actual  word  7777777777777776  (octal)  in  all  900 
words.  The  file  is  read  on  logical  unit  number  3.  E.ich  record  contains 
data  for  100  cells,  each  of  which  has  the  following  format: 

Word  1  Sector  number  .Format: 

LLLLLLLL  nimmnimmmm,  where 

L  -  Latitude  of  lower  right  corner  in  seconds,  fixed  point. 

m  -  Longitude  of  lower  right  corner  in  seconds,  fixed  point. 

2  State  and  city  code  (BCD)  Format: 

SSCCCCCC,  where 

S  -  State  code,  BCD 
C  -  City  code,  BCD 

3  City  center  and  SMSA  indicators  (BCD)  .Format: 

OOOOOOCS  where 

C  -  city  center  indicator 
S  -  SMSA  indicator 

4  Nighttime  population,  fixed  point 

5  Daytime  population,  fixed  point 

6  Industrial  output,  floating  point 

7  Area  of  sector,  floating  point 

8  Population  density  of  sector,  fixed  point 

9  Number  of  tracts  in  sector,  fixed  point 
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Each  cell  read  is  tested  to  see  if  the  population  density  is  greater 
than  MINDEN  anti  if  the  industrial  density  is  greater  than  MININD.  If 
either  is  above  the  cutoff  the  cell  is  stored  into  memory;  if  not,  the  cell 
is  accumulated  into  one  of  10  "tail  cells"  according  to  its  maximum  pop¬ 
ulation  density.  MINDEN  and  MININD  are  read  in  the  basic  data  deck. 

If  sense  switch  3  is  on,  the  program  gives  a  printout  of  the  stored  cell 
data  ir:  a  form  corresponding  to  the  above  format.  This  can  be  of  two  types, 
depending  on  the  variable  NAMEC,  the  last  card  in  the  basic  data  deck. 

If  NAMEC  is  blank ,  all  cells  will  be  printed.  If  NAMEC  is  some  city  code, 
then  only  that  city  will  be  printed. 

When  the  terminal  record  is  recognized,,  unit  3  is  rewound,  and 
control  returns  to  the  option  selector. 

Run  Type  LVCTDAT 

This  is  the  option  which  reads  in  the  populatior  tract  data  for  use  in 
the  targeting  model.  Again  the  data  file  is  a  single-reel,  binary  buffered 
file,  this  time  with  996  words  per  phyaical  record.  The  file  ir  unlabeled 
with  the  actual  word  7777777777777776  (octal)  in  the  terminal  record.  This 
file  is  read  on  logical  unit  number  4.  Each  record  contains  data  for  83 
tracts  in  the  following  format: 

Word  1  Sector  number.  Format: 

LLLLLLLLmmmmmmmm  where, 

L  =  Latitude  of  lower  right  corner  in  seconds,  fixed  point 
m  =  Longitude  of  lower  right  corner  in  seconds,  fixed  point 
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Word  2 


Tract  latitude  in  degrees,  floating  point 

3  Tract  longitude  in  degrees,  floating  point 

4  Standard  location  code  (BCD) 

5  City  center  and  SMSA  indicators  (BCD).  Format; 

OOOOOOCS  where 

C  -  city  center  indicator 
S  -  SMSA  indicator 

6  City  code  (BCD) 

7  1965  daytime  population,  fixed  point 

8  1965  nighttime  population,  fixed  point 

9  1970  nighttime  population,  fixed  point 

10  1970  daytime  population,  fixed  point 

11  Area  of  sector,  floating  point 

12  Industrial  output,  floatingpoint 

The  city  whose  code  is  NAMEC,  read  in  the  basic  data  deck,  will  be 
that  used  for  the  data  for  the  targeting  model. 

Run  Type  MINIMAX 

This  option  enables  one  to  control  whether  defenses  are  to  be  generated  using  the 
mixed  lambda  (or  single  lambda)  procedure  or  the  minimax  procedure.  Both  procedures  are 
defined  in  the  description  of  subroutine  DEFOPT.  Tie  minimax  approach  tries  to  limit 
the  excess  kill  over  what  could  be  achieved  for  each  of  several  lambdas. 

Only  one  parameter  needs  to  be  set  to  change  to  minimax  operation;  that  is  the 
switch,  MINIMAX  --  1.  It  is  suggested  that  a  smaller  number  of  mus  be  used  and  that 
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NNEW,  the  number  of  mus  used  on  second  and  subsequent  passes,  be  of  the  order  of  five 
or  ten  for  program  efficiency.  The  input  definitions  are  in  that  section,  and  sample  values 
for  these  are  shown  in  Sample  Run  Deck  *2. 

R  un  Type  OPTIONS 

This  run  type  is  used  to  set  one  of  the  parameters  IDPOP,  IAP0P,  IFILL,  or  IAOB 
to  any  value  desired.  For  the  format,  see  the  data  sequence  definition;  for  examples  see 
the  sample  run  decks. 

Run  Type  PERTURB 

This  option  changes  the  lethal  areas  according  to  any  desired  assumed  change  in 
hardnesses,  although  the  change  should  be  small,  since  partial  derivatives  are  used  to 
change  the  areas.  The  parameter  which  determines  the  change  is  DELTA,  the  positive  or 
negative  fractional  change  in  hardnesses.  After  reading  this  parameter,  the  lethal  areas 
are  adjusted,  and  a  revised  lethal  area  table  is  printed.  PERT  URB  ,is  performed  in  sub¬ 
routine  PERTURB  in  OVLAY3. 

Run  Type  TRIMTAIL 

This  option  enables  one  to  set  NCELLS  to  any  desired  value.  It  can  be  used 
for  such  things  as  evaluating  a  defense  without  the  tail  cells  being  included. 

3.  Output  Options 

Run  Type  OUTPUT 

This  option  enables  one  to  obtain  additional  information  about  (he  defense  and 
the  attack  on  the  defense.  The  routine  does  three  breakdowns— defenses  vs.  population 
density,  attack  vs,  density,  and  attack  vs.  defense.  Figure  llshows  how  the  population 
density  categories  are  set  up,  with  the  number  of  people  in  each  category.  Next,  in  the 
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same  figure,  an  important  breakdown  of  defenses  by  population  density.  It  is  interesting 
to  note,  for  example,  that  defense  options  3,  5,  7,  9,  14  and  15  were  never  chosen. 
Figure  12 shows  the  distribution  of  the  different  attack  densities  over  the  population 
density  brackets  and  over  the  defenses.  No  data  is  required;  the  option  card  alone  calls 
for  the  print.  This  routine  is  executed  in  subroutines  POPQUANT  and  CORELATE. 

Run  Type  PUNCH 

Inclusion  of  a  PUNCH  card  before  generating  a  defense  enables  one  to  obtain  a 
card  output  deck.  This  contains  all  the  defenses  assigned  over  the  country  for  future 
reference,  and  may  be  used  as  an  input  for  DEFSPEC.  No  data  is  required  for  this  option 
the  punching  is  executed  in  DEFPUNCH. 
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Figure  11.  Population  Density  Categories 
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4.  Auxiliary  Computation 

1 

1 

i 

Run  Type  CROSSCOR 

j 

This  option  sorts  the  cells  into  population  and  industry  density  brackets 


and  keeps  tallies  of  the  marginal  distributions  as  well  as  the  joint  distribution. 

The  parameters  to  set  are  indicated  in  the  data  definition  section.  The  population 
density  brackets  are  given  in  population  per  square  nautical  mile;  the  industry 
density  is  thousands  of  dollars  per  square  nautical  mile,  and  the  marginal  distribu¬ 
tions  for  these  brackets  are  included  in  the  heading.  The  correlation  table 
accumulates  industry  in  each  population  density- industry  density  "box"  and 


normalizes  by  the  total  amount  of  industry  in  each  industry  density  bracket. 


Run  Type  EXERCISE 


This  option  enables  one  to  obtain  tuned  (single  lambda)  defenses  for  a  number 
of  cost  levels  and  lambdas,  evaluating  each  of  them  for  all  attack  objectives.  More 
importantly,  it  then  uses  the  lambda  which  generated  the  defense  to  find  an  attack  against 
the  defense,  thus  obtaining  a  "solution  to  the  Lagrangian"  in  the  sense  of  the  theoretical 
section  above.  Otherwise  its  usefulness  is  mostly  for  small  scale  data  in  exploration  on 
the  behavior  of  different  lambdas  for  different  costs. 

Before  entering  EXERCISE,  all  the  same  parameters  need  to  be  set  as  in  SIMPLE; 
in  addition,  the  attack  objectives  must  be  set  (presumably  using  ATTOBJ).  EXERCISE 
uses  subroutine  WORKOUT,  which  reads  an  array  of  costs,  and  of  lambdas,  and  of 

attack  levels.  For  each  cost  and  defense  lambda,  DEFOPTG  is  called,  and  then 
DEFPUNCH,  as  well  as  EVAL  and  ANSPRIN  for  every  attack  objective.  Then  for  every 
attack  objective  the  attack  using  the  defense  lambda  is  obtained,  in  OFFOPTZ,  and 
that  "solution  of  the  Lagrangian"  is  printed. 

Run  Type  MiXLAM 

This  option  enables  one  to  observe  the  effect  of  changing  the  relative  weights  of 
two  lambdas  in  generating  defenses.  The  routine  explores  defenses  plus  evaluations  for 
combinations  of  the  two  lambdas  in  steps.  Its  use  is  in  small  data  analysis,  for  the 
computation  time  is  very  large  otherwise.  The  routine  first  tries  an  (unsuccessful)  method 
of  guessing  the  mixed  lambda  combination  and  then,  if  desired,  scans  all  possible  combi¬ 
nations  of  the  two  lambdas  in  steps  of  .05.  After  each  defense  is  found,  an  evaluation 
is  performed,  and  an  output  is  made  which  is  identical  to  SIMPLE.  For  input  data,  see 
the  input  sequence  definition  section. 
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Run  Type  TABLDENS 

This  option  tables  the  optimum  attack  densities  and  payoffs  as  for  different  values  of 
p/\  ,  where  p  is  the  population  (or  payoff)  density.  This  data  is  for  information  only — it 
is  not  used  in  the  program.  Figure  13shows  a  sample  table  for  the  first  defense  option. 

For  example,  for  f>/\  less  than  .312  but  greater  than  .102,  the  optimum  attack  density 
is  .06  and  the  payoff  is  .615  with  the  height  of  burst  being  6  (air  burst).  This  printout 
is  repeated  for  all  the  defense  options. 


.men  as-sims  ro<  tatr  ■  1 


attack  8E*stTy 
o.oo 

tor  RhO/LA^HA  LESS  THAN 
0.13167 

attack  payoff,  lower  density 
o.ooooi 

“Oft  attack  payoff,  upper  density 

0  .  _ 0.61471 

MCI 

-  6 

0,06 

0.13 

0.31234 

0.9*050 

0.61471 
....  0.01482 

6  0,61402 

6  6.94229 

6 

6 

0 ,25 

0,50 

5.Q2160 

37.15274 

0,94229 
_  0.99200 

6  0,99208 

0  0.99060 

4 

0,?5 

1,90 

2*6.67702 

2005.03 997 

0.99000 

0.99982 _ 

6  0,99902 

1  .  1.00000 
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6  1,00006 
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00 
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1.00000  : 
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1  -  _ U0OOO0  .. 

1 

Figure  13.  Sample  Table  for  the  First  Defense  Option 
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C.  GLOSSARY  OF  VARIABLES  IN  COMMON 


Basic 

A  (12) 

Value  (to  attacker)  of  each  payoff  component 

ACOST 

Actual  cost  spent  if  defense  saturates 

ALPHA 

Weight  on  offense  cost  in  defense  Lagrangian 

AREA 

Current  cell  area  (sq.  n.mi) 

AREAMAX 

Normalizing  constant  for  cell  area 

ATTCORR(20 ,  15) 

Frequency  count  for  15  defense  options  and  20  attack  densities 
in  CORELATE 

ATT  LEV  (20) 

Attack  levels  (in  equivalent  megatons) 

ATTPAY(IOO) 

Total  attacker  payoff  for  each  current  lambda 

ATTPAYZ(20) 

Best  estimate  of  attacker  payoff  for  each  attack  level 

ATTSPRED(10,20) 

Frequency  count  for  10  density  brackets  and  20  attack  densities 
in  CORELATE 

BASEPOP 

Base  population  selected  by  IDPOP. 

BDLAM 

Lambda  for  current  balanced  defense  calculation 

BDPSI 

Shelter  hardness  assigned  to  current  cell  for  balanced  defense 

CDEF  (300) 

Cost  of  optimum  defense  for  each  mu,  total  over  cells 

CDEFCL  (300) 

Cost  of  optimum  defense  for  each  mu,  current  cell 

CDLAG  (300) 

Optimum  Lagrangian  for  each  mu,  current  cell 

CELLCOST 

Balanced  defense  cost  for  current  cell 

CITPOP  (650) 

City  population 

COSTPER  (15) 

Cost  per  person  of  defense  option  (shelter  mix) 

Preceding  Page  Blank 
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CSTSHEL  (10) 


D  (20) 

DCOST 
DEFEXCL 
D LAM  (22) 

ERMF.AS  (20) 
EXCLUDEF(15 ,  15) 

FDEP  .10,15) 

FIND 

FINDMAX 

FKIL  (12,20,15) 

FLAM  (100) 
FLAMLEV  (20) 
FLAMX(IOO) 
FLAMZ(IOO) 

FMU  (200) 
FMUX(30G) 

FRACP  (300) 
HLAG(IOO) 

HOB  (15) 


Cost  per  shelter  space  of  10  shelter  types. 

Possible  attack  density  quanta  (equivalent  megatons/sq.n.m.) 
Defense  cost  in  dollars  for  current  defense  option 
=  1  if  current  defense  is  excluded  (set  by  DCONSTR) 

Lambdas  considered  by  the  defense 

Error  measure  (fraction  of  attack  level)  for  each  lambda  in  EVAL 

Used  in  time-phasing;  for  each  possible  final  defense  option 
EXCLUDEF  is  1  for  each  excluded  intermediate  defense  option, 
otherwise  0 . 

Fraction  of  population  supplied  one  of  10  shelter  types  for  each 
of  15  defense  options. 

Industrial  value  of  current  cell 

Normalizing  constant  for  cell  industrial  value. 

Fraction  of  12  hardness  components  killed  by  20  attack  densities 
at  15  heights  of  burst. 

Initial  set  uf  offense  multipliers  for  evaluation 
Best  estimate  of  lambda  for  each  attack  level 
Set  of  offense  multipliers  for  any  call  to  0FF0PT 
Temporary  array  for  new  lambda  generation 
Values  for  current  iteration 
Initial  mu  values 

Fraction  of  split  cells'  areas  on  primary  defense 
Optimum  offense  Lagrangian  for  each  lambda,  current  cell 
The  15  possible  heights  of  burst 
Optimum  payoff  for  each  lambda,  current  cell 


H0PT  (100) 


IAOB 

iAPOP 

ICBEG 

ICEND 

ICELL 

IC0MP(12) 

ICT1  (650) 

ICT  2(650) 

IDEF 

IOENS 

IDFOPT  (300) 

IDOPT(IOO) 

IDPOP 

IFiLL 

I H I  (20) 

I  HID  LAM 
IHOB 

IHOBOPT  (100) 
ILODLAM 
ILSAVE 
ISAT 

I5PLC  (300) 


Index  of  attack  objective 

Population  basis  at  attack  (second  and  third  population  bases): 

1  -  night,  2  -  day,  3  -  maximum,  4  -  minimum 

Index  of  first  cell  of  group  to  be  optimized 

Index  of  last  cell  of  oroup  to  be  optimized 

Current  cell  index 

Hardness  level  of  et:ch  (non-empty)  payoff  component 

Starting  index  in  cell  data  for  city 

Ending  index  in  cell  data  for  city 

Index  of  defense  shelter  option 

Attack  density  index 

Index  of  optimum  defense  for  each  mu,  current  cell 

Optimum  attack  density  index  for  each  lambda,  current  cel1 

Population  basis  for  defense  option  (first  population  basis): 

1  -  night,  2  -  day,  3  -  maximum,  4  -  minimum 

Shelter  filling  mode:  1  -  optimal,  2  -  partial  alert, 

4  -  balanced  defense 

Index  of  bracketing  lambda  (large  of  two  lambdas)  for  each  attarx 
level  in  EVAL 

Index  of  lowest  (variable  misnamed)  lambda 
Height  of  burst  index 

Optimum  height  of  burst  indes  for  each  lambda,  current  cell 
Eauivalent  to  IHIDLAM 

Index  of  closest  (lower  est.)  mu  in  cost  closing 
-  2  if  defense  cannot  spend  TCOST,  otherwise  =  1 
Cell  index  of  split  cells 
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1ST  1  (60) 


IST2  (60) 

ISW1 

ISW2 

JDEF 

JDEFSP  (300) 

JPOPDENS(IO) 

KDEF 

MAPHOB  (15) 

MINDEN 

MINIMAX 

MlNIND 

MMAOB 

NAHOB 

NALEV 

NAMCIT  (650) 
NAMSTAK60) 
NAOB 

NCELA(5000) 
NCELB  (5000) 
NCELLS 
NCIT 


Starting  index  in  cell  data  for  state 

Ending  index  in  cell  data  for  state 

Setting  both  =  1  allows  changing  cell  data 

Normalizing  constants  FINDMAX,  AREAMAX,  POPMAX 

Defense  option,  current  cell 

Secondary  defense  option  of  split  cells 

Population  density  brackets  in  POPQUANT 

Auxiliary  defense  option  (for  constraint  purposes) 

Indices  of  heights  of  burst  scanned 

Population  density  cutoff  for  cells  stored  in  memory 

~  1  for  minimax  selection  of  cell  defense  Lagrangian  over  lambdas 

Industry  density  cutoff  for  cells  stored  in  memory 

Originally  for  minimax  over  attack  objectives,  not  now  used 

Actual  number  of  heights  of  burst  scanned 

Number  of  attack  levels 

City  names  (BCD) 

State  names  (BCD) 

Number  of  attack  objectives 

Packed  ceil  data,  to  be  unpacked  by  RDCELL 

Number  of  population  cells 
Number  of  cities 
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NCOMP 


NDEF 

NDENS 

NO  LAM 

NHOB 

NLAM 

NLA  MX 

NMU 

NMUX 

NNEW 

NPOPDEN(IO) 

NPSI 

N5PLIT 

NSPLTOLD 

NSTAT 

NSTYP 

PAYIND(IOO) 

PAYINDZ(20) 

PAYPOP(lOO) 

PAYPOPZ(20) 

PKILL(IOO) 

PKP 


Number  of  hardness  components  in  payoff  (non-empty) 

Number  of  defense  options 

Number  of  attack  density  quanta 

Number  of  lambdas  considered  by  the  defense 

Number  of  height  of  fcurst  options 

Number  of  attacker  Lagrange  multipliers 

Number  of  current  offense  multipliers 

Number  of  initial  mu  values 

Number  of  mus  in  current  iteration 

Number  of  mus  used  in  second  and  subsequent  passes  in  DEFOPTG 
Number  of  cells  in  density  bracket  in  POPQUANT 
Number  of  hardness  levels 
Number  of  split  defense  cells 

Number  of  split  cells  in  previous  defense  (for  time-phasing) 

Number  of  states 
Number  of  shelter  types 

Total  industrial  value  destroyed  for  each  current  lambda 
Best  estimate  of  industrial  kill  for  each  attack  level 
Total  population  kill  for  each  current  lambda 
Best  estimate  of  population  kill  for  each  attack  level 
Optimum  population  kill  for  each  lambda,  current  cell 
Daytime  population  of  current  cell 
Maximum  population  of  current  cell 


PMAX 


POPACT 

POPCOUNT(IO) 

POPMAX 

PCPMULT(IO) 

POPSH(IO) 

POPUL(12) 

PSK15) 

PSIMAX 

PSIMIN 

PUNCH 

RES 

SHLSPRED(10 , 1 5) 

SPACES(IO) 
SPREAD  (12,10) 

SPRIND  (12) 

STATE 

STPGP(60) 

SWITCH 
TCOST 
TOLERR 
VALIND  (12) 


Actual  current  cell  population 

Population  in  density  bracket  in  POPQUANT 

Normalizing  constant  for  cell  population 

Factor  to  convert  population  kill  to  attacker  value 

Population  sheltered  in  10  shelter  types,  current  cell 

Population  in  each  hardness  component 

The  15  possible  ps  levels 

Upper  cutoff  psi  for  balanced  defense 

Lower  cutoff  psi  for  balanced  defense 

Settinq  -  1  gives  punch  output 

Nighttime  (resident)  population  of  current  cell 

Frequency  count  for  10  density  brackets  and  15  defense 
options  in  CORELATE 

Spaces  available  in  10  shelter  types,  current  cell 

Fraction  of  population  in  each  rtf  10  shelter  types  assured 
vulnerable  at  each  of  12  hardnesses. 

Fraction  of  industrial  value  at  each  of  12  hardnesses 

Setting  -  1  gives  STATCOST  defense 

State  population 

Variable  read  by  option  selector 

Target  cost  of  defense  investment  for  group  of  cells 

Tolerable  percent  error  in  attack  level 

Industrial  value  in  each  hardness  component 
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VINDMULT(IO) 

VKILL(IOO) 

W(15 ,15) 

WDAOB(lO) 

WDLAM(15) 

WEPEQ(lOO) 

WMIN 

WPEXP(IOO) 


Factor  to  convert  industrial  value  killed  to  attacker  value. 

Optimum  industrial  kill  for  each  lambda,  current  cell 

Lethal  area  for  each  height  of  burst  and  hardness 

Weight  defense  assigns  to  each  attack  objective 

Weight  defense  assigns  to  each  lambda 

Optimum  weapon  expenditure  for  each  lambda,  current  cell 

Lethal  area  corresponding  to  upper  psi  cutoff,  balanced  defense 

Total  weapon  expenditure  for  each  current  lambda 
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Glossary  (continued) 
Bounding  Procedure 


BCOST(3> 

BKILL<3) 

BLAM(3) 

3MAX 

BMIN 

BMU(3) 

BWEP(3) 

CUMC0SK20 ,20) 

CUMKILL(20  ,20) 

CUMWEP(20,20) 

CUR  COST  (20  ,20) 

CUR  LAG  (20 ,20) 

CURWEP  (20  ,20) 

EPS 

HILAM 

HSAT 

JJCH0Z(22) 

MODESW 

NENO 

SATWEP(22) 
SLOPE  (20) 


Cost  at  beginning,  middle,  and  end  of  iambda  Sweep 
Payoff  at  beginning,  middle,  and  end  of  lambda  sweep 
Lambda  at  beginning,  middle,  and  end  of  lambda  sweep 
Upper  bound  values 
Lower  bound  values 

Mu  at  beginning,  middle,  and  end  of  lambda  sweep 

Weapon  exp.  at  beginning,  middle,  and  end  of  lambda  sweep 

Cumulative  cost  for  20  lambdas  and  20  mus 

Cumulative  payoff  for  20  lambdas  and  20  mus 

Cumulative  weapon  expenditure  for  20  lambdas  and  20  mus 

Current  cell  cost  for  20  lambdas  and  20  mus 

Current  cell  Lagrangian  for  20  lambdas  and  20  mus 

Current  cell  weapon  expenditure  for  20  lambdas  and  20  mus 

Parameter  controlling  separation  between  lambdas  chosen 

Maximum  lambda  which  need  be  considered 

Maximum  payoff  possible 

Index  of  closest  mu  for  each  lambda 

Selects  mode  of  running  verification 

Number  of  end  limits  reached  on  lambda  sweep 

Weapon  expenditure  at  saturation,  line  upper  boundary  payoff 

Slope  of  line  lower  boundary  payoff 
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VERIFY 


Setting  -  1  turns  on  verification  procedure 
YCEPTLO  (20)  Intercept  of  line  lower  boundary  payoff 

YNTCEPT  (22)  Intercept  of  line  upper  boundary  payoff 

D.  INPUT  DATA  SEQUENCE  DEFINITIONS 

Tli i 3  section  contains  the  definition  of  the  data  sequences  for  the  various  run 
options  or  types  and  for  the  basic  data  deck;  the  data  must,  of  course,  appear  in  this 
form .  The  data  for  each  option  begins  with  the  option  card,  shown  at  the  top  of  each  de¬ 
finition,  and  continues  or  branches  according  to  the  defined  sequences. 

The  formats  are  all  exactly  as  they  appear  in  the  various  subroutines;  the  subscripts 
on  all  subscripted  variables  are  dummy  variables  in  the  program  and  are  shown  here  as 
I  or  J . 

In  addition  to  the  options  given  here,  options  EVAL,  LIVEDATA,  LVCTDAT 
OUTPUT,  PUNCH,  STOP,  TABLDENS,  TIMEFAZE,  and  UNSHEL  have  no  data 
cards  except  the  option  card  itself.  Option  STATCOST  and  SIMPLE  have  the  same 
data  sequence,  but  STATCOST  allocates  money  to  states  by  population  and  then 
optimizes  the  nationwide  defense  state  by  state.  Option  WIGGLE  has  been  dismantled 
and  if  used,  has  the  same  effect  as  a  SIMPLE  card. 
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BASIC  DATA  PECK  INPUT  SEQUENCE 


FORTRAN  Variable  Format 

NOENS  (8U0) 

D<0  (8F10.2) 

NHOB,  NPSI,  LACOMP  (8110) 


HOB(J)  (8F10.2) 

PSI(I)  (8F10.2) 

If  LACOMP  p  0  skip  to  *  in  sequence . 

WCJ.I)  (8F10.2) 

Skip  to  **  in  sequence. 

W(J,I)  (8F10.2) 

NABOB  (8110) 

MAPHOB(J)  (8110) 

NSTYP ,  NOEF  (2,10, 

FDEF  (I ,J)  (8F10.2) 

SPREAD!!  ,J)  (8F10.2) 

SPRINO(I)  (8F10.2) 

CSTSHEL(I)  (8F10.2) 

FINOMAX,  AREAMAX,  POPMAX  (3F10.2) 


Explanation 

Number  ol  attack  density  quanta. 

Density  quantum  levels. 

NHOB  is  number  of  heights  of  burst  (  A  15) . 
NPSI  is  the  number  of  hardnesses  (  4  15). 

If  LACOMP  p  0,  lethal  areas  are  computed  from 
lethal  radii . 

Here  read  NHOB  heights  of  burst. 

Here  read  NPSI  hardnesses. 


Here  read  lethal  areas  (sq.  n.mi)  by  HOB  (on 
single  cards)  and  PSI . 


Here  read  lethal  radii  (kilofeet)  by  PSI  (on  single 
card)  and  HOB. 

NAHOB  is  the  actual  number  of  heights  of  burst 
to  be  scanned  (4  NHOB). 

If  NAHOB  =  NHOB  skip  next  card  in  sequence. 

MAPHOB(J)  are  the  indices  of  the  NAHOB  hoglit 
of  burst  scanned. 

NSTYP  is  the  number  of  shelter  types. 

NDEF  is  the  number  of  defense  options. 

Here  read  fraction  of  each  defense  option  J  composed 
of  each  shelter  type  I  by  defense  option  (on  single 
card)  and  shelter  type. 

Here  read  fraction  of  each  shelter  type  J  at  each 
hardness  I  by  shelter  type  (on  single  card)  and 
hardness. 

Fraction  of  Industrial  Value  assumes  to  have  each 
of  NPSI  hardnesses. 

Cost  per  person  for  each  ol  NSTYP  shelters. 

Maximum  industrial  value,  area,  and  population 
over  all  cells . 


ICELL,  RES,  PKP,  FIND,  AREA, 

JDEF ,  KOEF  (110 , 4F10 .2  ,  2110)  ICELL  is  the  cell  number,  RES  is  nighttime  popula¬ 

tion,  PKP  is  daytime  population,  FIND  is  industrial 
value,  AREA  is  area,  JDEF  is  defense  option  and 
KDEF  is  for  constraint  purposes.  Continue  reading 
these  cards  (cell  data)  until  ICELLSO. 


NALEV 
ATTLEV(I) 
MINOEN ,  MININD 


NAMEC 


(8110)  Number  of  attack  levels  (  <  20). 

(8F10.2)  Here  read  NALEV  attack  levels. 

'2110)  Lower  population  density  and  industrial  density  cut¬ 

offs  for  a  cell's  inclusion  into  memory.  A  cell  is 
included  if  MAX  |  PKP.RESl/  AREA iMINDEN  or 
FIND/AREA  ^MININD. 


If  NAMEC  is  blank,  the  cell  printout  in  option 
LIVEDATA  will  list  all  cells;  if  it  is  a  city  code, 
the  listing  will  be  of  that  city  only. 


Fortran  Variables 


Dau  SwjuefKt  fof  Run  Type 


^ ALERT 


Format 


EypiaivUtort 


falert 


rar  io.2) 


Fr^ct.'on  of  population  aietted . 


Data  Sequence  for  Run  Type  (  ATTOBJ 
I 


Fortran  Variables 

Fomat 

£*nlanaiion 

NAOB 

(110) 

Number  of  attack  objectives 

POPMULTlIl 

(8F10.2) 

Here  read  population  kill  fnult«plie*s 
for  NAOB  attack  objectives  . 

VINDMULTtll 

(8F10.2) 

Industrial  kill  multipliers. 

WOAOBII) 

(8F10.2) 

Here  read  weight  defense  assigns  to 

NAOB  attack  objectives. 

Data  Sequence  for  Run  Type  |  BALDEF 

Fotuan  Variable 

Format 

» 

Explanation 

OEFNAMt. NEWOAIA 

(A8,2X,I10) 

DEFNAME  is  any  8  character  code  for  the  balanced  defense 

If  NEWDATA  =  0  skip  to  ***  m  sequence. 

The  following  lethal  area  data  (down  to  w')  needs  to  be 
specified  for  unsheltered  posture  only. 

NHOB,  NPSI,  LACOMP 

(8110) 

NHOB  is  number  of  heights  of  burst  (  ~  15). 

NPSI  is  the  number  of  hearnesses  ( •“  15), 

If  LACOMP  p  0,  lethal  areas  are  computed  from  lethal  radi 

HOB  (J) 

(8F10.2) 

Here  read  NHOB  heights  of  burst. 

PSItl) 

(8F10.2) 

Here  read  NPSI  hardnesses. 

II  LACOMP  t  0 

skip  to  *  in  sequence. 

WtJ,l> 

Skip  to  **  in  sequence 

(8F10.2) 

Here  read  lethal  areas  (sq.  nm  )  bv  HOB  ton  single  cards) 
and  PSI. 

WU,I) 

(8F10.2) 

Here  read  iethal  radii  (kilofett)  by  PSI  (on  single  card)  and 
HOB. 

NAHOB 

(8110) 

NAHOB  is  the  actual  number  of  heights  of  burst  to  be 
scanned  (  -  NHOB). 

IF  NAHOB  -  NHOB  skip  next  card  in  sequence. 

MAP  HOBIJ) 

(8110) 

MAP  HOBtJ)  are  the  indices  of  the  NAHOB  heights  of  burst 

SCAflMd. 


BALOEF,  Cont. 

Fortran  Variables 

Format 

Explanat  ion 

NSTVP  =  2 

i 

(8110) 

Two  shelter  types  —  unsheltered  and  variable  hardness  sheltered 

NDEF  =  2 

(8110) 

Two  defense  options. 

SHLNAME,  CSTSHEL,  NHARD,  XK 

(A8,2X,F10.2, 

I10,F10.2) 

Here  read  unsheltered  dati 

SHLNAME  is  the  name  of  the  unsheltered  posture. 

CSTSHEL  is  the  cost  pet  person  of  the  unsheltered 
posture  (usually  =  0.0) 

NHARD  is  the  number  of  hardnesses  of  the  unsheltered 
posture. 

If  NHARD  =  1  then  XK  is  the  hardness  of  the  unsheltered 
posture  in  psi  and  skip  next  card. 

If  NHARD  ^  1  then  XK  may  be  blank  and  read  next  card 
to  obtain  hardness  spread.  Hardness  must  appear  in  PSI 
array. 

HARD  (1),  FRAC(l),  HARD  (2),  FRAC  (2), . . .  (8F10 .2)  Harnesses  of  unsheltered  posture  in  psi  and  fraction  of 

population  assumed  to  have  each  hardness. 

SHLNAME,  CSTSHEL,  NHARD,  XK  (A8,2X,F10. 2,  Here  read  variable  hardness  data. 

I10,F10.2)  SHLNAME  is  arbitrary. 

CSTSHEL  is  arbitrary  or  blank  (it  is  not  used). 

NHARD  =  1. 

XK  should  be  one  of  the  hardnesses  used  for  the  unsheltered 
posture. 


I  1 

it — r~ 

1  l 

(3110) 

These  two  cards  define  the  defense  options  for 
balanced  defense. 

2j 

l|  2| 

»  ! 

»  i 

- 

(3110) 

1 

SPRIND(I) 

:  ! 

(8F10.2) 

Fraction  of  industrial  value  assumed  to  have  each  of 

NPSI  hardnesses. 

NCOST 

(110) 

Number  of  costs  for  which  defenses  will  be  computed. 

COST  (1) 

(-9PF10.2) 

Here  read  NCOST  costs  (in  billions  of  dollars). 

NALEV 

(110) 

Number  of  attack  levels. 

ATTLEV  (1) 

(8F10.2) 

Here  read  NALEV  attack  levels  (in  equivalent  megatons). 

PSIMIN,  PSIMAX,  WMIN 

(8F10.2) 

PSIMIN  is  the  lower  PSI  cutoff.  If  balanced  defense  assigns 
hardness  below  this  levei,  then  unsheltered  posture  is  assigned 

PSIMAX  is  the  maximum  allowable  hardness  for  a  shelter. 
WMIN  is  the  lethal  area  corresponding  to  PSIMAX. 
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Data  Sequence  fof  Run  Type 

^  BOUND 

f"  j  r:  ^ 

£»p‘dn^*Of> 

MODESVV 

T10J 

Th»$  selects  a  '-ode  for  the  bounriioc 
procedure.  After  initial  Jala  go  to 

Sequence  for  node  desired. 

NEWVij 

(110. 

ff  NEWMU  *  0  read  new  »’••«{..  ntherwis? 
skip  to  rode  sequence  desired. 

fmumj 

■8  F10.2; 

Here  read  20  new  rous. 

/'  3 

Mudt,  | 

1 

This  node  computes  initial  point  from-  a 
cost  and  a  lambda,  then  sweeps  lambdas  fror 
10,0  to  1000000.0  to  generate  bounds. 

DLAMi  1  i 

(F10.2) 

This  is  the  initial  Iwnbda. 

TCOST 

(F20.2) 

Total  cost  of  shelter  posture  (input  irt  doll^s). 

HSAT 

(F20.2) 

Total  payoff  (or  population). 

EPS 

(F10.2) 

Parameter  controlling  separation  between 
lambdas  chosen.  Approximately  it  isthe 
distance  between  upper  and  lower  bounds. 

Made  j'  2 

This  mode  reads  in  initial  point,  otherwise 
as  m  Mode  1 . 

BKILL;2),  SWEP(2>,  BOOST !2), 
BLAMI21,  BMU2) 

(F20  2,  F10  2,  F20  2, 
2F10.2) 

Initial  kill,  weapon  expenditure,  cost, 
lambda,  and  nu. 

HSAT 

(F20.2) 

EPS 

(F10.2) 

/~\ - - 

Mode  j 

This  mode  reads  in  initial  point;  then 
lambdas  are  chosen  to  verify  this  in.ua>  point, 
and  bounds  are  Set  for  the  attack  lewei  of  the 
initial  point.  The  data  sequence  is  exactly  as 
in  Mode  2 , 

Mode  |  ^ 

This  mode  computes  initial  point  from  a  cost 
and  an  attack  level ,  then  verifies  that  point 
as  m  Mode  3  . 

ATTLEV  (11 

(F10.2) 

Attack  level  for  initial  point. 

TCOST 

!F20.2) 

HSAT 

(F20.2! 

EPS 

(F10.2) 

Mode  f7” 

This  mode  reads  lambda s  and  does  one 
pass,  computing  the  bounds  for  those 
lambdas. 

TCOST 

(F20.2) 

NOLAM 

(110) 

Number  of  le  ■  *.  l  &  20i 

OLAMfl) 

(8F10.2I 

Here  reai  Wl‘L  AM  lambdas  ‘ascending  order) 

HSAT 

!F20.2) 
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Mode 


Mode 


V 


*-  orrrrn  vr 


t*p'inar 


BXILL12),  BWEP(2),  BC0ST(2), 
8LAMI21,  BMU(2) 

T  COST 

NDLAM 

DLAM(I) 

HSAT 


(F20.2,  F10  2.F20.2, 
2F10.2! 

(F20.2) 

(110) 

(8F10.2) 

(F20.2) 


BKILUl),  BWEP(l),  BCOST(l),  4(016, 4X) 

BLAM(l),  BMUIl ) 

BKILL<3),  BWEP(3),  6C0ST(3),  4(016,4X) 

BLAM(3),  BMU(J) 

TC0.'>T,  HSAT,  EPS  (-9PF  10.2,  OPF  20.2, 

-6PF10.2) 


This  "'Ohe  rp.irts  0  point  t-i  hi*  ver.lp.  1  thor. 
w"  pules  I.w..jk!s  (or  thji  v, ; l  aud  ,  !,  v,_.; 
using  lambdas  which  are  read  m. 

Kill,  weapon  expenditure,  cost,  l.ipihu.. . 
and  mu. 


Number  ol  lambdas  (  *  20) 

Here  reud  NDLAM  lambdas  (ascending  order) 


This  mode  lor  use  when  sweep  type  run 
(1  of  2>  is  terminated.  Uses  octal  dump 
to  continue  computation. 

Lower  attack  level  end  ol  previous  sweep. 


Upper  attack  level  end  of  previous  sweep. 


Fortran  Variables 


Data  Sequence  for  Run  Type 


^ CR0SSC0R 


Format 


Explanation 


NBOXES,  IDPOP,  INORM,  (81101 

1CBEG,  ICENO,  MIN 


NBOXES  is  the  number  of  correlation 
categories  (S 31). 

IDPOP  is  the  population  basis  =  1,2. 
or  3  lor  NIGHT,  DAY,  and  MAX  . 

INORM  is  the  normalization  mode 

=  1  nonnalized  by  cells  in  population  bracket 
-?  normalized  by  cells  in  industry  bucket 
2  normalized  by  population  in  bracket 
=4  nonnalized  by  industry  in  bracket 

ICBEG  is  tne  beginning  index  of  the  quup  of 
cells;  ICEND  is  the  end. 

MIN  is  the  lower  cutoff  ol  density  brackets; 
all  cells  below  tins  density  are  placed  in  one 
category. 
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Lr?  •.♦r,  --.e  •*  *  ?  .-•» 


r  sc-f  ' 


*■  ortra-  var<rb*e 


'0" 


KEY,  TCOST,  NCEILS 


( 10X  ,110 . 20K  , -9PF  1  0  2  ,  KEY  IS  toy  integer  nu  - Mr  to  idenl.tr  the  de‘ense 
10X.HO)  TCOST  it  the  cost  of  the  defense 

NCELLS  is  the  nu- aer  of  celis  included  I*  NCELLS  r 
old  NCELLS  on  error  message  is  printed  end  rrma.ninq  cells 
are  set  to  unsheltered  condition 


!DSP£C  (!> 


NSPUT 


(40125 


(110) 


Here  read  NCELLS  defense  options,  one  lor  each  cell  at  40 
per  cond.  EachlOSPEC  £  NDEF  except  for  split  ceils, 
where  index  is  100  -  defense  option, 

Nu-ber  o‘  solit  cells 


ISPLC(I),  JDEFSP(I).  FPACP(I)  (2110, F10  8)  Here  mad  NSPLIT  split  cell  data  cards —  index  ol  split 

cell,  secondary  defense  ol  spilt  cell,  and  fraction  on  primary 
defense. 


If  this  card  is  GO  continue  in  sequence. 

If  this  card  is  SKIP  ,  skip  to  *  in  sequence. 
If  this  card  is  QUIT  skip  to  ***  m  sequence. 


IDPOP,  IAP0P,  IFILL,  ISAT,  (6110, 2L10) 

NDLAM,  NAOB,  MINIMAX.  MMAOB 


IDPOP  and  IAPOP  are  population  bases,  IFILL  is  shelter 
filling  mode,  ISAT  is  saturation  indicator,  NDLAM  is  number 
ol  lambdas  lor  defense,  NAOB  is  number  ol  attack  objectives, 
arid  MINIMAX  and  MMAOB  deal  with  minimax. 


TCOST, ACOST, ALPHA 


DLAM(I) 


(-9PF  10.7,  10,7, 
0PF10.8) 


(8F10.2) 


WDLAM(I),  WDA0B(1),P0PMULT(1), 
VINDMULTtl  ),WDA0B(2l,  P0PMULTI2) 


(8F10.8) 


TCOST  is  the  shelter  cost 

ACOST  is  actual  cost  in  case  of  saturation 

ALPHA  Is  normally  1,0. 

Here  read  NDLAM  lambda  considered  by  defense. 

Here  read  NDLAM  weights  for  lambdas  and  NAOB  attack 
objective  weights,  population  multipliers,  and  industry 
multipliers. 


NHOB,  NPSI,  LACOMP  <8  1 10) 


HOB(J) 

PSI  (0 

II  LACOMP  ^  0  skip  to  +  in  sequence. 


(8F10.2! 

(8F10.2) 


If  this  card  is  GO  continue  in  sequence. 

If  thi 5  card  is  SKIP  skin  to  *■’-  in  sequence. 

If  this  card  is  QUIT  skip  to  *+*  in  sequence 

kfctethat  IA0B  is  not.  an  input  variable  in  this 
data  sequence.  It  may  be  put  in  using  the 
OPTIONS  sequence 

NHOB  is  number  of  heights  of  burst  (¥15). 

NPSI  is  the  numbet  of  hardnesses  (  5  15). 

II  LACOMP  X  0,  lethal  areas  are  computed  from  lethal 
radii. 

Here  read  NHOB  heights  of  burst . 

H*»e  read  NPSI  hardnesses. 


57 


Ot*  <Pf  : 


Portr*t!  v'  »ia*?r 

r  t~r 

f  if  {H  tnm  pftrt 

W(J,I> 

Skip  to  -t*  -+-  in  sequence  . 

8F10  2) 

Mere  read  lethal  area;  'sq.  nr  )  by  ”08  ion  tingle 
car ds  and  PSi 

+-  W(J,I) 

(8F 10 .2  i 

Here  read  ieLhil  radii  (kilofeet)  by  PSI  (an  single 
card)  and  HOB. 

-i- NAHOB 

(8  1  10 

NAHOB  ts  the  actual  number  ol  hearts  a)  bury 
to  be  scanned  (  <  NHOB). 

II  NAHOB  -  NHOB  skip  next  card  in  sequence. 

MAPHOB(J) 

(8  1  10) 

MAPHOB(J)  me  the  indices  of  the  NAHOB  height 
ol  bust  scanned. 

/go”  /new  /skTp  /quTt 

If  this  card  is  GO  continue  in  sequence. 

1  wj  °M  "1 

If  this  card  is  NEW,  insert  entire  INPUT  sequence 
(without  INPUT  card).  Note  that  with  NEWDATA  '  0 
it  isnot  necessary  to  repeat  RDATT  sequence  No  more  data 
required  after  that . 

If  this  card  is  SKIP  or  QUIT  skip  to  ***  in  sequence. 

*»  NSTYP ,  NDEF 

<21101 

NSTYP  is  the  number  of  shelter  types. 

NDEF  is  the  number  of  defense  options. 

FDEF  (l,J) 

(8F10.2) 

Here  read  fraction  of  each  defense  optiun  J  composed  of  each 
shelter  type  1  by  defense  option  (on  single  card)  and  shelter 
type. 

SPREAD  <I,J> 

(8F10.2) 

Here  read  fraction  of  each  shelter  type  J  at  each  hardness 

1  by  shelter  type  (on  Single  card)  and  hardness. 

SPRIND  (1) 

CBF10.2) 

Fraction  of  Industrial  value  assumed  to  have  each  of 

NPSI  hardnesses. 

CSTSHEL  (1) 

(8F10.2) 

Cost  per  person  for  each  of  NSTYP  shelters. 

/noeval  /  ”” 

If  this  card  is  NOEVAL,  no  evaluation  will 
be  performed,  otherwise  DEFSPEC  will  evaluate 

1  o»  1 

the  defense  with  parameters  as  set . 

n  c  «  n  T  ["EXERCISE 

Data  Sequence  for  Run  Type 

Fortran  Variable 

Format 

Explanation 

NCOST 

0X01 

Number  of  costs  far  which  defenses  will  be  computed. 

COST(I) 

(4F20.2) 

Have  read  NCOST  costs  (in  dollars). 

NWLAM 

(110) 

Number  of  lambdas  for  which  defenses  will  be  computed. 

WEAMtl) 

(8F10.2) 

Have  read  NWLAM  lambdas. 

NALEV 

(110) 

Number  of  attack  levels. 

ATT  LEV!  II 

(BF10.2) 

Have  read  NALEV  attack  levels  (in  equivalent  megatons). 
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L.  l  t  *V»  '  -  i  ^i'4i*AcA 


OEFNAMc.  is  any  h  character  tone  io»  shf‘l«  input 


■*o,  M,  4  0i 


NhuB,  NPi  ,  LACOMP 


'6  I  10) 


HOB  <J) 
PSI  <11 

If 

A’!,  1,1) 


LACOMP  A 


0  skip  to  * 


18F10.2) 

<PF10  2> 
in  sequence. 

(8F10.2> 


Skip  to  r*  in  sequence. 


II  NCWDATA  r  0  skip  to  *"*"*  in  sequence. 

NHOB  is  number  ol  heigh's  ol  burst  <  4  1 5 » . 

NPSl  is  tbe  number  of  hjrpnc-sses  '  4  15). 

If  LACOMP  X  0,  lethal  areas  me  computed  from  lethal 
radii . 

Here  read  NHOB  heights  of  burst . 

Here  reed  NPSI  hardnesses  . 

Here  read  *eth.£i  areas  (sq.nmi  by  HOB  (on  single  cards) 
and  PSI. 


'  WU,I> 
NAHOB 


(8F10.2I  Here  read  lethal  radii  tkilofeet)  by  PSI  (or  sing):  card)  and 

HOB. 

(8110)  NAHOB  is  the  actual  number  of  heights  of  burst  to  be 

scanned  (  4  NHOB). 

If  NAH0B=  NHOB  skip  nest  card  in  sequence. 


MAPHOB  (J) 


(8  110) 


MaPHOB  (J)  are  the  indices  of  (he  NAHOB  height  of  burst 
scanned. 


N5TYP 

of 

these 

groups 


NSTYP 

NDEF 

"SHLNAWE,  CSTSHEL,  NHARD,  XK 


(  I  10) 

(  I  10) 

(A8,2X,F10,2, 

I10.F10.2) 


HARD  (1),  FRAC  (1),  HARD  (2),  FRACI2), . . .  (8F10.2) 


Number  of  shelter  types.  I  i  10) 

Number  of  defense  options.  (’'IS' 

SHLNAME  is  the  name  of  the  shelter-  any  8  characters. 
CSTSHEL  is  the  cost  per  person  of  the  shelter 
NHARD  is  the  number  of  hardnesses  of  the  shelter. 

If  NHARD  =  1  then  XK  is  the  hardness  of  the  shelter  in  psi 
and  skip  next  card. 

If  NHARD  1 1  XK  may  be  blank  and  read  neid:  card  to  obtain 
Hardness  spread.  Hardness  must  appear  in  PSI  array 

Hardnesses  of  shelter  in  psi  and  fraction  of  shelter  assumed  to 
have  that  hardness. 


NDEF 

of 

these 

groups 


IDEF,  NSHEL,  KX  (3  1  10)  IDEF  is  the  index  of  the  defense  option. 

NSHEL  is  the  number  of  shelters  included  in  the  option. 

If  NSHEL  -  1  then  KX  is  the  index  of  the  shelter  type  and 
skip  next  card.  If  NSHEL  r- 1  KX  may  be  b!a.-i«  and  read  next 
card  for  sheiter  mix. 

ISHEL  (1),  FRaC(1),ISHEL(2),FRAC(2),  . .  .(-1(110, FI 0.2)>  Indices  of  shelter,  types  and  fraction  ef  mix  composed  of  each. 


SPRIND  (I) 


<8F  10.2)  Fraction  of  industrial  value  assumed  to  have  each  of 

NPS i  hardnesses. 
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D  -St -H  fV  T  .  : 


FarUau  Vtii*uu>g 

MINIM  AX,  MMAOB.  NDLM, 
NMLO 


FMU  (I) 


F  U*-1  «M  £p£^4Mj<KV> 

MINIMAX  1  «.!l  cause  DEFOPT  to  mr.  m  ---  m.  • 
{2L10.6I10V  procedures;  0  lor  normal  procedure. 

MMAOB  for  minima*  attack  objectives,  not  used  *  - 
should  he  blank  or  zero 

!f  NDLM  0  10  -  ore  data  is  rsQuiicd.  If  NDLM  >s  «or- 

zero,  it  becomes  NMU  and  ne*t  card  is  read 

If  NMLD  is  nonzero  it  beco^s  NNEW,  the  ru-ber  of  i-n's 

cn  t^e  second  vast  »n  DEFOPTG. 

(8E10.4)  Here  read  NM U  mus 


Data  Sequence  for  Run  T  ype  ^ MIXLAM 

Fortran  Variable  Format  Explanation 


XMLAM(l),  XMLAM(2)  (2F10.2)  Two  lambdas  to  be  used  for  mixed- 

lambda  combinations. 

TCOST  (-9PF10.2)  Total  cost  of  defense  posture  (in 

millions). 


If  this  card  is  EXPLORE,  program  will 
explore  mixed- lambda  combinations  m 
stops  of  0 .05  after  trying  first  mixed- 
lambda  guess. 


I''' EXPLORE  or 


Fortran  Variables 


Data  Sequence  lor  Run  Type 


'''options' 


Format 


Explanation 


SWITCH,  (VALUE 


r 

r 

r 

r 


IDPOP 


IAP0P 

or 

IF  ILL 
or 


(A8 , 2X ,  110) 


SWITCH  is  examined  to  see  if  it  is 
IDPOP,  IAP0P,  IFILL,  or  IA0B. 

If  it  is,  the  program  sets  whatever  variable  «l 
is  to  IVALUE  and  reads  the  next  card.  If  it 
is  not,  die  run  type  terminates,  and  SWITCH 
is  used  as  the  next  run  type  card. 


Dali  Sequence  for  Run  Type  ^  PERTURB 

Fortran  Variables  Format  Explanation 

DELTA  (8F10.2)  Fractional  change  in  hardnesses  of 

sheitm  (positive  or  negative) . 
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>  a Iran  Vqiabit 


ICOST 

NNDLAM 


OLAM(I) 

NWDLAM 


.  t’  i 


>  t  O-.J#  f 


t  W-M 

(•WiO.JI 

(110) 


(8F10.2) 

(110) 


EiplaoJtiQn 

Total  cou  of  shelter  pasture  finput  'n  billions' 

Number  of  lambdas  to  V  considered  b y  de'ense  (  «  15) 

l>  NNDLAM  -  0  prolan  u$ej  previous  lambdas;  skip  to 
*  sequence 

Here  read  NNDlAM  defense  lambdas. 

If  NWDLAM  -■  0  program  assigns  equal  weights  to  lanbdas; 
skip  to  »  in  sequence. 


WDt.  AM 

(6F10.2i 

NAOB 

1110) 

POPMULT(I) 

(8F10.2 

VINDMULTfl) 

(8F10  2) 

I  NOEVAL  or 

1 

NEWDaOB 

1110) 

WDAOBl  1) 

(8F10.2) 

have  read  NNDLAM  weights  for  defense  iwnbdas. 

Number  ol  attack  objectives  f  5  10).  If  NA0B  =  1, 
no  more  data  >s  iwimwt 

Here  read  popu  at'on  kri  rnultipliers  for  NAOB  attack 
objectives 

NACB  industrial  kill  multipliers.  After  reading  this 
card  the  program  generates  the  defense . 


If  NOEVAL  card  appears,  no  more  date  is  reguired  and 
this  run  type  terminates  w  (haul  evaluating  the  defense. 

If  NEWDAOB  -  0,  no  more  data  is  required.  Otherwise, 
read  ne*t  card. 

Here  read  weights  defense  assigns  to  NAOB  *tack 
objectives.  Program  generates  new  defense  and 
evaluates  it  for  all  attack  objectives. 


Return  to  **  in  sevence . 


Fortran  Variable 


Data  Sequence  for  Sun  Type 


Format 


j/TRIMTAIL 

Explanation 


NCELLS 


NCELLS  is  the  number  of  cells  considered  in  the  run. 
The  number  may  be  cut  down  to  exclude  tail  cells,  or 
may  be  set  for  any  otter  purpose. 
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E.  VALUES  OF  IMPORTANT  PROGRAM  CONSTANTS 


The  following  variable,  are  initially  set  in  the  program  to  the  values  indicated. 

These  may  then  be  reset  as  desired  using  the  input  sequences. 

IDPOP  -  l 

IAPOP  =  1 

IFILL  =  1 

IAOB  ~  1 

NAOE  =  1 

NDLAM  =  1 

POPMULT(I)  =1.0 

VINDMULT(I)  =  0.0 

WDAOB(I)  =  1.0 

DLAM(I)  =  2500.0 

WDLAM(I)  =  1.0 

TCOST  =  10000000000.0 
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III.  SAMPLE  DECK 


A.  BASIC  DATA  DECK 

Most  of  the  basic  data  deck  exists  only  because  the  program  grew  over 
a  period  of  time,  and  this  method  of  reading  data  once  and  then  writing  over  it 
insured  that  runs  did  not  terminate  because  a  parameter  was  unset.  There  is  no 
reason  to  ever  change  cards  on  most  of  this  data  deck;  the  important  variables 
(the  ones  which  are  not  normally  overwritten)  will  be  noted  so  changes  may  be 
made. 
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First  is  read  the  number  of  attack  density  levels  and  the  density  quantity 

5 

t 

themselves.  These  density  levels  (especially  the  minimum  nonzero  one)  affect  j 

the  way  the  density  model  approximates  the  damage  effects;  the  values  here  were 

used  on  all  the  FINAL  shelter  runs  except  one.  This  data  is  not  overwritten,  ' 

A 

and  the  attack  density  levels  are  set  only  in  the  basic  data  deck.  \ 

i 

§ 


I 


These  cards  indicate  one  height  of  burst  (  =  1234.5)  and  10  psi  levels,  and  give 
the  quantities.  The  hardness  data  and  small  population  data  in  this  basic  data  deck  correspond 
to  the  data  used  in  the  analytic  model  in  Appendix  C. 
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The  10  lethal  areas  for  the  simple  height  of  burst  are  next  (these  were  computed 
from  L.A.  =  121/psi). 
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This  is  the  actual  number  of  heights  of  burst. 
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Here  is  indicated  7  shelter  types  and  7  defense  options. 
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This  is  the  FDEF  array;  there  are  7  defense  options  consisting  of  pure 
deployments  of  the  7  shelter  types. 
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This  is  the  SPREAD  array,  indicating  7  pure  shelter  types  at  3,  10, 
30,  70,  100,  200,  and  300  psi. 
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This  is  the  SPRING  array  -  the  definition  of  industrial  hardness  for  this  basic 
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Here  are  the  shelter  costs  and  then  FINDMAX,  AREAMAX,  and  POPMAX 
for  the  small  cell  data  which  follows. 
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Here  is  the  12  cell  data  used  in  Appendix  C  ;  this  time  with  area  in 


nautical  miles.  There  are  twelve  cells  and  for  each  cell  is  read  ICELL,  RES,  PKP, 
FIND,  AREA,  JDEF,  and  KDEF.  The  twelve  are  terminated  by  the  blank  card. 

This  data  is  useful  for  checkout  and  exploratory  runs. 


69 


STAMCARD  i  WORD  v&*»D 


All  the  previous  data  except  the  attack  density  levels  is  normally  reset  to  desired 
lethal  area,  shelter,  and  cell  data.  The  following  cards  are  not  normally  reset  and  serve 
as  parameters  for  all  programs  run. 
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Here  are  tfie  12  attack  levels  (in  equivalent  megatons)  used  for  evaluation  and 


verification  runs. 
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These  are  parameters  MINDEN  and  MININD,  the  population  density  and  industry 


density  cutoffs.  When  the  full  scale  data  is  read,  a  cell  is  included  explicitly  (rather 
than  in  the  tail  cells)  if  either  its  population  density  is  greater  than  7,000  per  sq  .mi. 
or  its  industry  density  is  greater  than  3,500  per  sq.mi. 
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This  card  terminates  the  basic  data  deck .  Its  use  is  in  the  cell  listing  for  option 
LIVEDATA  —  if  the  card  is  blank  all  cells  vv.ii  Ue  listed,  if  the  card  contains  a  city 
code  (such  as  22530005  for  Washington,  D.C.)  only  the  cells  in  that  city  will  be 
listed. 


71 


QMI  •  OnvQtiViC 


B.  SHELTER  INPUT  DECK 


This  is  the  deck  which  determines  the  shelter  parameters  to  be  used  in  the  run. 

The  deck  is  optional  in  that  if  this  entire  deck  is  removed,  the  program  will  run  properly  using 
whatever  shelters  have  been  input  with  the  basic  data  deck.  Thus,  it  is  possible  to  simply 
use  the  basic  data  deck  to  set  the  shelter  data.  However,  the  particular  form  of  input  used 
in  that  portion  of  the  program  is  inconvenient,  especially  if  changes  in  shelter  mixes  or 
definitions  are  made.  All  runs  in  the  study  were  made  using  this  INPUT  option,  rather  than 
the  basic  data  deck,  to  put  in  shelter  data. 

The  shelter  input  deck  shown  here  is  the  data  comprising  the  FINAL-2  shelters,  which 
were  used  for  the  bulk  of  the  study.  For  a  complete  description  of  what  is  allowed  for  the 
INPUT  option,  see  the  section  which  defines  the  data  sequences  for  the  different  options. 
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This  card  causes  the  option  selector  to  give  control  of  the  program  to  the  INPUT 
option  portion  of  the  program.  After  going  through  the  entire  data  sequence  for  the  INPUT 
option,  control  will  again  return  to  the  option  selector. 
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On  'his  card  FINAL-2  is  a  name  for  this  shelter  input  deck.  The  1  is  a  switch 
which  here  indicates  that  new  hardness  data  will  be  read. 
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Th  >  card  indicates  10  heights  of  burst  and  12  hardness  level.  The  1  is  a  switch 
indie  tes  that  lethal  radii  will  be  read  in  kilofeet  with  all  heights  of  burst  for  a  single 
hardness  c.:  a  single  pair  of  cards. 
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These  are  the  12  hardness  levels,  ranging  from  5.4  psi  to  490  psi. 


fcc  .  V 

j  c  1  .  <r 

VsXRD  l 

nj  w  iro  2  g 

WORLi  3  f) 

WORD  *  H 

WORD  5  g 

WORD  6  g 

I’jObd  .  Til  w  ,Ho  a 

ni 


v  or:j  3  n  ord 


\  ord 


i 


Hi  H| 

il  o  o  on  o  o  oooo  005  0  0  0  ol  S  o  a  o  i  o  i  ojo  o:o  o  o  oo  o  o  ojb  b  o  oToToTol 

. • . — - . - 1 . -  -  -  .i*** 

I  t  l! 


i'l  .  i  ill  .  I  nmit.»»»»;iiii»i»!'rn»r. s'xaitajH 

. I.  ......  <l'i  .  .  i  i'i  i  i  i!  i'i  i  i  t  i  1  I  ii  pi  I  I  M 


I  I  1  7  1M  7  I  I 
7.7  2  2  ll?  2  7  2 
3!!  3  3  313  33 


1  ';i  1 1 1.1 1 nj'n  1  ill  11 1  l'.i  luim 

[2  7!?  2  2  2|2  2  2  2I2  2!?  7  2  7.2  2  2  jil  2  ?  2  7  21 7  7  2 

ll  J'I  3313  3  3  7 

1  ' 

4  4  4  <4  4'4  <  4  4 


3  31  3  3'3  3  3  3  3  3'l  3  3  31)  3  3  3 

I  1  '  I  *  1 

!<  4  4  414  4  4  4  4  4'4  4  4  4.4  4  4  4  «  l!<  4  4  <  4  4  4  4 

*  '  I  '  I  »  '  I  ,  I  I  I 

5  5  5  Sis  5  5  5  5  5  5  5  1  5  5  5  5  5  5  5'5  5  5  5|5  5  5  SI  5.5  5  5  51  5  5  5,5  5.5  5  5  51 


!??1?  7  2  2  ? 
13;  3  3  31 

1  I  ' 

4  II  4  4  4 ' 4 


lit  list t I 


t  t  j  1 1  ill  t  i  i  t  t  t  ititti  i  t  t  e  1 1 1>  1 1 1 


7  ,1  1  7  11  77  7  1  uniml?  7  77.7 


tititt't 


?|7  7  7  7,7  7  7  7|7  7;7  7  7  7|7  7  7  7 

111  1  I'l  I  I  tit  II  .  I  I  1  1  I  I  I  I  I'l  111  I  I  ill  l;:*  I  I  I  I  1 1  II  l!  •  I  I  III 

,1,  ..J.,,  9,9999  99  9  9  9  9  9  9  t!t  9  5  5  9  5  9  5  5  31  91  9  9  51  5  9  5)9  919  I  111  1 1  ill  1  9  9  9l9  9  9  9,9  9  5  9  5  91  9  S3. 

i  Jin  *11  J  ’a'a:.3*'n  uw  . »  »  »  j:  *  1.*.  .../«-«•  .1 ««  »ar  v  u»»  . .  ’  *  "*> 


V  CBD  O  gl  V  ODD  7 

ol 


,w.t . . toFoloTillii 

j»i  wu  v  »  tt*»»  1a  tt  *•••»  u|u  »•  •  tt|*J  ei  ••  tt  ti  n,n  h  n  »iif  •  tt  ■ 


to  o;o  0  0  oro  0  0  Vo  00  0  000  0  0  0 


ill  1 1  i!i  1 1 

31'!  3  3  3'3  3  3  Sj 

4  4*4  4  4  4  4  4  4 

■  . 

5  515  5  5  5:5  5  5 

Stitt  IS  II 

ill  7  ill  7  1 


2  ill??  217  7  7  2 


i .1  n  t.i  1 1  i  i  :i  1 1  r.t  1 1 1 


t  tit  i  ill  1 1  i'i  i'i  iilm  ill  i‘"i  i  tin  i1  I  J 

I  '  ■  i  .  1  i  '  '  r 


71  7  m  7  1 7 


3  3;  3  3  313  3  3  3  3  Jl  3  nil  J  3 1| 

4  4  4  4  4  4'4  4  4  4  «  4<4  4  4  4>4  4  4  4 

,  1  |  I 

5  5l5  5  51  5  5  5  5  51  5  515  5  5  5 

‘  »  I  1 

llllf  tl.llll  11,11111111 

I  1  1 

7  ?1?  7  7  ?J7  7  7  7  7  7J7  T  7  7,7  7  7  7 


These  are  the  lethal  radii  for  a  one-megaton  weapon  for  5.4  psi  hardness 
and  for  the  10  different  heights  of  burst.  The  lethal  radii  for  the  other  11  hard¬ 
nesses  follow  in  the  same  pattern,  10  heights  of  burst  for  a  single  hardness  on  a 
pair  of  cards,  then  moving  to  the  next  hardness.  These  are  shown  in  Table  I,  and 
they  follow  immediately  after  the  preceding  two  cards. 
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Table  I  .  Remainder  of  Lethal  Areas 
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Here  it  is  specified  that  only  two  of  the  10  heights  of  burst  will  be  scanned 
for  attack  optimization  -  numbers  1  and  6,  or  0  and  5500  feet. 
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Here  are  indicated  5  different  shelters  (  National  Fallout  Shelter  Survey, 
10  psi  blast,  30  psi  blast,  100  psi  blast,  and  300  psi  blast)  and  15  different 
defense  options  or  shelter  mixes. 
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Here  is  the  definition  of  the  hardness  spread  for  the  NFSS  shelter.  The 
name  of  the  shelter  is  NFSS,  the  cost  is  zero,  and  the  number  of  hardnesses  is  3. 
The  second  card  specifies  the  3  hardnesses  and  the  fraction  of  the  shelter  assumed 
to  have  each  hardness;  here  one-third  each  at  5.4  psi,  10.0  psi,  and  at  19.5  psi. 
The  hardness  specifications  for  the  remaining  four  shelters  follow.  Note  that  for 
each  of  these  shelters  10  psi,  etc. ,  is  just  a  name  for  the  shelter,  and  that  the 
characteristics  of  the  shelters  are  specified  by  the  hardness  spread  and  thence  by 
the  lethal  radii  previously  read. 
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Here  die  shelter  mix  specifications  begin.  The  first  defense  option  is 
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The  first  number  1  is  the  index  of  the  defense  option.  The  second  number 
is  the  number  of  she1  ter s  in  the  mix  for  that  defense  option.  If  there  is  only  one 
shelter  in  the  mix,  then  the  third  number  specifies  which  of  the  shelters  just  read 


is  the  one  shelter  type.  Here  for  the  first  defense  option  there  is  one  shelter  in 
the  mix  and  that  shelter  is  the  NFSS  type,  shelter  number  1. 
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For  the  second  defense  option  there  are  ?  shelter  types  in  the  mix.  When¬ 
ever  the  number  of  shelters  in  the  mix  is  more  than  one,  a  second  card  is  read 
which  specifies  the  mix.  Here  the  mix  is  0.8  of  shelter  type  1  (NFSS)  and 
0.2  of  shelter  type  4  (100  psi).  The  shelter  mixes  for  the  other  13  defense 
options  are  specified  similarly,  as  ihown  in  Table  II;  they  follow  immediately 
after  the  above  two  cards. 
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Table  II.  Remainder  of  Shelter  Mixes 
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The  last  two  cards  define  the  hardness  spread  for  industry  by  specifying 
the  fraction  of  the  industry  which  has  each  of  the  12  hardness  levels.  Here 
industry  is  assumed  to  have  the  same  hardness  spread  as  the  NFSS  shelter  type, 
or  one-third  at  each  of  5.4  psi,  10  psi,  and  19.5  psi.  These  two  cards 
terminate  the  data  sequence  for  this  INPUT  option.  Control  now  returns  to  the 
option  selector  for  whatever  type  of  run  is  desired. 
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C,  SAMPLE  RUN  DECK  #1 

The  purpose  of  this  run  is  to  generate  a  defense  using  the  mixed  lambda 
technique  and  assuming  a  mixed  attack  objective  and  then  to  evaluate  the  delens-, 
twice  for  two  different  attack  objectives.  The  defense  in  this  sample  deck 
found  with  IDPOP  =  1  and  IAPOP  =  2  with  the  budget  level  at  10.0  bill: 

The  following  data  sequence  is  put  behind  the  basic  data  deck  and 
shelter  input  deck: 
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This  card  causes  the  program  to  read  the  population  data  tape.  Control 


is  then  returned  to  the  option  selector 
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These  two  cards  set  IAPOP  =  2  (the  program  sets  IDPOP  =  I). 

The  first  card  sends  control  to  the  portion  of  the  driver  which  examines  the  sub¬ 
sequent  card  for  one  of  the  parameters  IDPOP,  IAPOP,  IAOB,  and  IFILL.  If 
the  card  is  one  of  those  parameters,  then  the  parameter  is  set  to  the  value  in 
column  20  and  the  next  card  is  examined.  If  it  is  in  the  list,  the  value  is  set; 
otherwise  it  is  assumed  to  be  the  next  card  for  the  option  selector  part  of  the 
program . 
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This  card  sets  a  switch  which,  after  the  defense  is  generated,  causes  the 
defense  options  selected  to  be  punched  onto  cards  for  every  cell.  Also,  the 
parameters  and  shelter  data  for  the  defense  are  punched,  this  punched  deck  is 
then  a  suitable  input  for  further  evaluation  runs  as  in  Sample  Deck  #  3.  After 
this  card,  control  returns  to  the  option  selector. 
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This  card  makes  the  program  do  run  type  SIMPLE,  which  is  the  basic  run 
type  to  generate  a  defense,  and,  if  desired,  to  evaluate  it.  The  entire  flexibility 
of  the  SIMPLE  run  type  is  described  in  the  data  sequence  definition  section;  the 
sequence  presented  here  is  one  possible  type  which  has  considerable  usefulness. 
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the  lambdas,  one  lambda  per  word  on  the  card. 
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These  two  cards  set  the  weights  for  the  lambda  just  read.  The  first  card 
may  be  any  non-zero  integer  (if  it  is  zero  the  program  weights  all  lambdas  equally 
and  the  second  card  will  not  appear).  The  second  card  contains  the  weights  for 
the  lambdas,  the  first  for  the  first  lambda  read,  etc.  The  mixed  lambda  combina¬ 
tion  just  read  is  78850 .0  weighted  at  0.6  and  7885.0  weighted  at  0.4. 
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Tliese  three  cards  read  the  attack  objectives  used  in  this  run.  The  ‘  irst 
card  is  the  number  of  attack  objectives,  NAOB.  If  NAOB  =  1,  then  the  run  is 
performed  using  the  attack  objective  set  previously  (the  program  sets  this  to  pure 
population  in  the  beginning).  Here,  two  attack  objectives  are  read.  The  first 
card  gives  the  population  multipliers  for  the  two  attack  objectives;  the  second 
card  gives  the  industry  multipliers.  The  first  attack  objective  is  a  50/50  attack 
objective  (equal  weightings  of  population  and  industry);  the  second  is  pure 
population.  The  first  weight  on  industry  is  1.077  and  not  0.5  because  there  is 
2.154  times  more  total  population  than  total  industry;  hence,  that  weight  is 
0.5  x  2.154  =  1.077. 
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STANCARD  8  tA GRO  CAR! 


When  the  program  generates  the  defense  it  will  assume  some  attack  objective 
for  the  attacker .  Since  the  weights  on  the  attack  objectives  were  set  to  1.0 , 0.0, 
0 .0 ,  0 .0 ,  . . . ,  the  defense  places  all  its  weight  on  the  first  attack  objective, 
the  50/50  mix.  Hence  this  defense  is  generated  assuming  a  50/50  attack  objec¬ 
tive.  The  second  attack  objective  is  not  used  in  the  generation  of  this  defense. 

At  this  point  in  the  data  sequence  Lite  program  will  compute  the  defense  and 
punch  it  onto  cards  for  future  use.  The  next  card  determines  whether  the  program 
will  immediately  evaluate  the  defense: 
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If  this  card  is  anything  but  the  key  word  N0EVAL,  then  the  program  will 
immediately  evaluate  the  defense  it  has  generated  for  all  attack  objectives,  one 
at  a  time.  There  it  evaluates  for  two  attack  objectives,  the  first  50/50,  and  the 
second  pure  population. 
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This  card  is  a  switch  as  to  whether  ihe  program  is  to  read  a  different  set  of 
attack  objective  weights  for  the  defense.  Here  the  zero  answers  negatively,  and 
til  is  card  terminates  the  SIMPLE  sequence.  If  the  card  were  any  other  number 
then  the  program  would  read  new  attack  objective  weights,  generate  the  defense 
and  evaluate  it.  Control  now  returns  to  the  option  selector. 
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This  card  is  the  normal  termination  for  the  program,  taking  the  program  out 
of  the  option  selector  loop. 
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Figure  14  shows  the  entire  data  sequence  for  this  run  type,  with  the 


basic  data  deck  and  the  shelter  input  deck  indicated  at  the  front. 


Figure  14.  Entire  Data  Sequence  for  Sample  Run  Deck  #1 
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D.  SAMPLE  RUN  DECK  #2 


The  purpose  of  this  run  is  to  generate  a  defense  using  the  minimax  tech¬ 
nique  and  assuming  a  mixed  attack  objective  and  then  tc  evaluate  the  defense  for 
both  day  and  night  pure  population  attacks.  The  defense  in  this  sample  deck  is 
found  with  IDPOP  =  1  and  IAP0P  =  2  with  the  budget  level  at  10.0  billion. 
The  following  sequence  follows  the  basic  data  deck  and  shelter  Input 

deck. 
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As  before,  the  first  card  causes  the  population  data  tape  to  be  read.  The 
second  card  causes  the  third  card  to  be  examined,  and  then  IAP0P  is  set  to  2 
(remember  IDPOP  =  1  already). 
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The  MINIMAX  card  here  enters  that  option  for  setting  up  the  minimax 
procedure.  The  second  card  sets  the  switch  MINIMAX  to  1  so  the  defenses  will 
be  chosen  according  to  that  procedure.  The  next  number  means  that  16  new  mus 
will  be  read  on  the  following  two  car  The  last  number  sets  NNEW  =  5; 

NNEW  is  the  number  of  mus  used  after  the  first  pass  in  DEFOPTG  —  the  value  5 
helps  speed  the  computation  better  than  leaving  it  at  50  as  it  would  otherwise  be 
set.  The  16  mus  read  next  are  tailored  to  fit  the  10  billion  budget;  they  need 
not  be  selected  thus,  but  it  will  help  computation  time  to  do  so. 
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Again,  this  card  causes  the  defense  generated  to  be  punched  out  for  saving. 


The  SIMPLE  option  follows  the  same  data  sequence  as  always  and  here 
sets  TCOST  =  10  billion,  and  reads  the  two  lambdas  to  be  used. 
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This  card  causes  WDLAM  to  be  set  to  0.5  for  each  lambda  (equal  weight). 
Even  though  the  lambda  weights  will  not  be  used,  the  data  sequence  must  be  main¬ 
tained,  and  this  is  the  simplest  method  here. 
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The  same  three  attack  objective  cards  as  for  Run  Deck  #1. 

At  this  point  the  program  will  generate  the  defense  using  the  minimax 
procedure  and  assuming  a  50/50  attack  objective.  It  will  then  be  punched  onto 
cards  for  future  use. 
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This  NOEVAL  card  terminates  the  SIMPLE  option  without  doing  the 
evaluation.  Control  returns  to  the  option  selector. 
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These  three  cards  cause  the  attack  objective  to  be  set  to  IA0B=2,  which 
is  pure  population  and  an  evaluation  to  be  performed  on  the  defense  just  gener¬ 
ated  for  the  pure  population  objective.  Now  the  point  of  jumping  out  of  SIMPLE 
and  doing  the  evaluation  this  way  is  simply  that  we  avoid  evaluating  the  defense 
for  the  50/50  attack  objective,  in  which  we  may  not  be  interested. 
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J  2  2  2  2  7  2  2  2  2-2  7  7  2'2  7,2  2  7  2-2  7  7  ?!:  2,2  2  2  7  7  7  7  7, 7  7  7  2  7  2  2  2  2  2i2  2^2  2  2  2;2  7  7  71  * 

mi’ i'i 1 1 i;i n ijni 1 1 1 1 n ill 1 1 333333  33 33 333133 1'3 1!]  1 3 33 3 3  3 


1  *  1  . 

4  4  4  4;4  4  4  4  4  4  4  4  4  4  4  4.4  4  4  4  4  4  4  1  4  4  I  4  4  4  4 

I'll''! 

5  5  5  5*5  5  5  5  5  5  5  5  5  5  3  5!5  5  5  5  5  3  5  5  5  5  5  5  5  5,5  5  5  5 


14444444444444444444 

1  t  j  i  i 

5  3  5  5  5  5  5  5  3  53  5'5  5  5  5l5  5  5  5 

6  S  5  6  5  6 B  S  i  6  B  S  5  G j  i  6  E  5  E  5  5  5  5  5  3  6  5  6  t  i  E  5  S'S  6  5  G  6  G  S  S  E  «'(  .  8  (  E  8  6  E  I  8 

7  7  7  ill  7;7  7  7  :  7  l  I  17  1  7  7  7  7  7  I  I  »|7  1.7  7  7  1,7  7  )  7  7  7  7  ’  7  7  7  7  7  l| 7  7;7  7  7  7,7  7  7  7j 


I  B  I  1  I  1  8  8  I  1  1  I  8  I  1  II  t  1  It  11 1 1  I  9  I  I  I  9  9  8  I  9  8  M  8  1 1 11 1 1  11  1 1  9'E  II  1  E  1 1  3  I  E9  1  3  I  9  3  I  1  I  I  I  HI  I  II  1  1 1 

1  i"1  9  9  9  9  9  9  9  9  9  9  1  9  9  9  9  9  9  9  9  9  9  9  3  1  9  9  9  9  1  3  9  9  9  9  9  9  9  9  3  9  9  9  9  9  1  9  1  ’  9  9  9  9  9  9  9  9  3  1  9  9  9  9  9  9  9  919  '3  9  9  9  9  9  9  9  9! 


■  ■;  u  u n  i»  1 1  is  n  .*«•?»  .ikss*;:  .« ®  r- 


i  i»  j*  k  *7  n  w  ei- i 


r.  U  «  .  IE  >M|  4«  «  V'  41  i;  »  *  IT  J*  SI  IS  U  IT  It  il 


1'9'J  9  9  9  9  9  9  9  91  I  J 

»-•  "4  *4  a  >  ’» I  r 


These  cards  change  the  time  of  attack  to  night  (IAPOP  =  2)  and  cause  a  new 
evaluation.  The  attack  objective  remains  at  pure  population,  and  all  other  parameters 
stay  fixed. 
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WORD  •» 


u 


III  WORD  r>  rjl  WORD 


WORD  O 


Ul 


DI 

9 , ,  3 , .  n  i  o  oToYooFTiioTa  oToTo'o  o'o  Vftfi  Sod  i  STitTiTSTTS  3  o‘o  o  if o  ITi  o'o  o  n  o  o  o  o  o  o  o  o  o'o  o  oTo  o’o  ST  o' 

*  •  i  i  ,  i  ■  *  !  id  i 1  "in '! 'i  ,  t  u  t* ::  jum n'r  .. n k  ’■  w*j  .  r  .4  ■  ■  .  «u>  «.mj “oki'hi'v  w  s  *mi  -.1  j;*'  cj  n  «  h  ti  w h  mih  a  *  n.f  *i  ^ i* 


;  II  I  I  1  1  1  1  1  I  1  I  1  It  1  I  III  M  1  1  11  1  1  1.1  111  1  1  M  1  I  I'l  Mill 


2  2  2  7  2  7’7  2  2  2;2  2!2  2  2  2'2  2  2  2 
.1  3  .  3  3  .  3  3  3  3  3  3  3  3  3  3)1  3  1  3 
4  ,  4  4  ,  4  4  4  4  4  4  4|l  4  4  ill  4  4  4 


2  2\2  2  2  2  2  2  2  212  2  2  2  2  2  2  2  7  2  2  7,2  2  2  2 

3  3  1  3  3  3'3  3  3  3l3  3  3  3  3  3  1  3  3  ill  3  3  3  3  3 


4  4,444  44444 


4444444444 


4  4  4  4  4  4 


5  5  3  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5)5  5'5  5  5  5  5  5  5  5  5  5i5  5  5  5  5  5  5  5:5  5  5555 

E  i  S  E  I  6  I  G  C  G  t;s  I  (  (!s  1 1 1!(  EE  G  S  BE  6  b  GIG  G'G  G  G  G  G  6  6  EE  BEEBE 

7  7  7  7  7  7  7  7!7  Tl  7  7  7  7  7  7  7,1  hi  111111111111  11  1  I  1|7  7-7  7  7  1 


1  I  I  l  I  11  1  1  I  1  1  1  I  1  I  1  1  I  U  II  I'l  1;l  I  1  I  :t  1  I  I 

2  2  7  712  2'2  2  2  7  7  7  7  2  2  2,2  2  2  2  7  2  2  2!2  2'7  2  2  2'2  2  2  l\ 

1  ■  i  ■  ■  i  1  ; 

3  3  n;i  3  3  3  3  3  3  3  3  3!3  33  3  3  3!l  1 1 3  3  l|l  1  1  3  3  3  llj 

4  1  1  1  4  4  4  4  4  4l4  4  4  ill  4  4  4  4  4'4  4  4  414  4'4  4  4  4  4  4  4  4!" 

;  '  " 

5  5  5  5:5  5,'5  5  5  5  5  5  5  5!5  5  5  5  5  5  5  5  5  5!5  5'5  5  5  5  5  5  5  5 

6  E  8  e!b  S'5  E  5  C  E  B  E  6  E  5  5  B  B  5  3  6  G  t  B  B  S  I  B  E  6  S  E  El 

'  '  :i  |  ; 

7  l  l  l|7  7!l  1  7  17  I  ?  7  7  7  1  1  7  7  7  7  ’  7  7  7  7  7  I  7  7  7  7  J!  s| 


i ;  9 1 4 1  a  1 1 1  a  1 1 1 1  it  1 1 1 1  n  1 1  m  1 1 1 1 1 1  a  1 1 H 1 1 1 1 1 1  a  m  i  is  1 1 1 1  i.i  1 1 !  s  t  *  3  •  9  ii  1 9 :  a  1 1 1 1 1 1 1  it 

9  3  9  9  9  9  9  9  9  9  3  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  919  1:9  9  9  9  9  9  9  9!9  9  9  9  a  9  9  9  9  919  9:3  9  9  9  9  9  9  9  9  9  3  5  5  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9 
4  3*119  -9.r  •?  1)  i»  is  u  •:  i:  s  x  :r  n  ”  :o  x  :t  :t  .1  n  1*  j)  jo  b  h  jr  ®  «:«  •««*«•»  u  m  y  a  *  j.  z  j  »i;ii  *:  u  h  4  **•  •  •;  :» r. 


The  OUTPUT  card  obtains  the  printout  described  in  the  OUTPUT  section, 
and,  finally  the  STOP  card  terminates  the  run. 
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E.  SAMPLE  RUN  DECK  #3 


The  purpose  of  this  run  is  to  read  in  a  previously  generated  defense  and  to 
evaluate  it  for  a  50%  alerted  nighttime  population  in  a  pure  population  attack.  The 
defense  to  be  evaluated  is  contained  on  a  punch  deck,  exactly  as  it  came  out  from 
some  previous  run.  Since  shelters  are  stored  on  that  output  deck,  we  need  have  only  the 
basic  data  deck  in  front  of  the  following  cards: 


WORD  2  [J  WORD  3  []  WORD  4  J]>  WORD  5  f)|  WORD  G  H  WORD  7  f|  WORD  R 


ill  !  iS 

"VI  •  ri  WORD  2  nj  WORD  3  [}•  WORD  4  WORD  5  [J  WORD  G  [j  WORD  7  WORD  Q  '  !| 

I  ‘  l  i 
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9  9  0  9  0  9  3  C  O  d'O  0  0  OTffTTO  O  OlToToTO  OOTOO  O'OOTO  0  0  0  9  0  0  0  0  0  0  0  oTTo' O'O  0  0  o'o  00  0  0  (Til  0  0  0  0  0  0  o  n  0  0  o'  ‘ 
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iiiiv  f  11 1 1 1 1  tit  11 1 1 11 1 1 1 11  nit  1:1 1  ii  !  ii  i|i  n  1 1  i;i  1 1  n  1:1  n  i;i  u  in  1111111111111111111  < 

j  1 2  j :  ?  j  1: 2 1 : 2 2 :  1  n  n;?  n  2  2  2 : 2!:  37  n  :•?  ?  ?  n  ii  j  j  17  1  ?  rln;?  2  j  2  2  ?  1  r?  2!* :  2  2  2  2  2  ii  n  2  o :  2  2  2  21 ! 

•  1  I  *  1  i  ;  ■  •  !  »  .  I  ?  j 

3  2  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3*3  3  3  3  3  3  3  3  3  3>3  3  3  3  3  3  3  3  3  3  3  3‘3  3  3  3  3  3  3  3l3  3;3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  2  3  3  3  3  3  3  3  3  3  3; 

4  1  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  «  1  4'-1'  «  4-4  4  4  4  4  4  4  4  4  :|l  4  4  4  4  l!l  4  1  4  4  4  4  4  4  4'l  4  4  4  4  4  4  4  4  4  4  1  4  4!  j 
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I  .  I  '  I  1 
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The  LIVEDATA  card  brings  in  the  population  data,  and  the  DEFSPEC  card 
causes  the  program  to  go  to  the  portion  which  reads  the  punch  deck  next: 


For  a  description  of  what  is  contained  in  the  punch  deck,  see  the  data  sequence 
definition  for  DEFSPEC. 
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Lit  woe,u  3  WORD  A  f  jj  WORD  5  [J  WORO  fc  f||  WORD  7  j]|  WORD 
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4  4  4  4  4  4  4  4  4  4  I  4  4  4.4  4.4  4  4  4j4  4  4  I  4  4  4  4  4  4  4  4  4  l|4  4.4  4  4  4.4  4  4  4j4  4 
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6  6  6  6  6  6,6  6  6  6.6  6  6  6|6  5  6  6  6  6  6  6  6  6|6  6  6  6  6  6  ,i  6  6  '  6  6  6  5  5  6 


5  5  6  6  6  6 


3  7  3  1:1 1,1  I  i  :  1  1  1  V,ll 1  1  1  VI  I  I  ljl  11  >  t  I  7  4  1  I,;  •  I  I  ;  )  7  )  I  1  1  7 


4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  '' 

'  1  ■  1  | 

5  5  5  5  5  5  5  5:5  5  5  5  5  5.5  5  5  5; 

I  !  1 

6  6  6  5  6  6  6  6  6  6.6  6  6  6  6  6  6  6  . 

7  I  I  1  1  7  7  7  7  7  7  7  7  7  7  7  7  7!  i 


I  1  1  5  1  9  8  1  1  8  5  4  3 1 1 I  I 1 1 I  I  l  a  8  1  I  S  1 1  t'l  I  I 1 1  4  1 1 1  lit  1 1 1 1 1 6  6  1  4  5  «  4  2  11  1  1 3  9  1  8  1  1  3  4  3  1  6  1  1  1  1  1  3  1  1  1  1  1; 
9  3  3  9  3  9  9  9  9  6  9  9,9  J  3  9  9  9  3  916  5  5  3  9  5  9  9  9  9!9  9:9  9  9  9  J  9  9  919  9  9  9  S  '  9  4  9  919  9  9  9  9  1  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  J  9  9  9  9  9, 
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This  NOEVAL  card  causes  DEFSPEC  to  not  evaluate  the  defense; 
otherwise,  the  defense  would  be  evaluated  with  parameters  as  they  were  set. 
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1 1 1 1 1  n  1 1  lit  11 1 1 1 1 1 1  i[i  11 1 11 1 1 1  iji  1:1 1 1  H  1 11:1 11 1 1 1 1 1  1  il  11 1 1 1,1  1 1  1 1 1 1 1 1  1 ;  1 1 1 1 1 1 1 1 11 1  1  |I° 
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These  cards  set  the  time  of  attack  to  night  and,  we  will  suppose,  the 
attack  objective  to  pure  population  (if  the  defense  we  are  evaluating  is  like  the 
previous  sample  decks,  this  would  be  pure  population).  Although  the  punch  deck 
saves  most  parameters  and,  indeed,  saves  the  POPMULT  and  VINDMULT 
arrays  (and,  hence,  the  attack  objectives),  it  does  not  save  the  particular  attack 
objective,  IAOB,  leaving  that  to  be  set  before  the  evaluation. 
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These  two  cards  set  up  the  fractional  filling  mode  and  define  the  fraction 
alerted  -  0.5.  All  parameters  are  now  set  for  the  evaluation. 
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Asa  last  step  we  evaluate  the  defense  using  the  new  parameters  and  conclude 


the  run. 
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F.  SAMPLE  RUN  DECK  #4 


The  purpose  of  this  run  is  to  find  the  upper  and  lower  bounds  for  the  optimal  curve 
for  $20  billion  defenses.  Following  the  basic  data  deck  and  shelter  input  decks  are  the 
following  cards: 
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These  cards  call  in  the  population  data  tape  and  set  IAPOP  =  2  . 
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These  cards  turn  control  over  to  the  bounding  routine  and  set  the  mode  switch  to  1. 
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These  cards  read  in  the  new  mus.  Notice  that  although  the  entire  range 
is  six  orders  of  magnitude,  a  number  of  the  mus  are  concentrated  in  what  is  the 
expected  region  for  the  $20  billion  level. 
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Here  the  initial  lambda  is  30,000,  TCOST  is  $20  billion,  and  HSAT 
is  209,473,668,  the  total  daytime  population  (since  no  attack  objective  was 
set,  it  remains  at  the  program  setting  which  is  pure  population  attack;  hence,  the 
total  payoff  is  the  total  population). 
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This  is  EPS ,  a  parameter  which  controls  separation  between  the  lambdas  in  the 
sweep,  and  hence  run  time.  This  value  of  2,000,000  was  used  in  all  runs  and  gave 
good  results,  although  some  runs  (as  this  one  did)  required  extra  computation  in  Mode  1 
for  desired  accuracy. 
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These  are  the  20  new  lambdas.  The  two  ends  (2034.69  and  30000.0)  were 
lambdas  used  in  the  computation  of  the  previous  case.  The  object  here  is  to  fill  in  a 
number  of  lambdas  over  the  same  region  as  before  to  improve  the  accuracy.  The  other 
lambdas  are  simply  spaced  between  the  ends. 

And  to  complete  the  sequence: 
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IV.  SUBROUTINE  DESCRIPTIONS 


The  program  is  divided  into  a  monitor  program  (contained  in  program  DRIVE 
and  the  overlay  programs)  plus  subroutines.  The  discussion  of  these  routines  is 
contained  in  this  section. 

One  brief  word  about  the  flow  charts  will  be  included  here.  Iterative  loops 
will  be  indicated  as  follows: 

V  Cel!  L _ _ 


The  symbol  V  means  "for  every,"  and  the  example  means  that  the 
indicated  loop  is  to  be  performed  until  all  cells  have  been  exhausted.  The  following 
indicates  a  table  lookup  process  of  locating  the  variable  FIND  in  the  table  RANGE: 
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f^°2  ram  DRIVE 

DRIVE  reads  the  basic  data  deck  (in  doing  so  it  calls  RDATT,  DEFINP,  and 
RDCELDT),  and  then  reads  option  cards  and  turns  control  over  to  the  proper  overlay. 
In  this  capacity  it  acts  as  part  of  the  option  selector;  the  same  function  is  performed 
in  the  first  overlay. 

Subroutine  ANSPRIN 

This  is  the  basic  answer  printing  routine;  its  output  is  discussed  under 
Run  Type  SIMPLE. 

Subroutine  BALDEF2 

This  is  the  driver  for  balanced  defense  computation.  It  is  similar  to 
DEFOPTG,  in  that  its  task  is  to  close  to  the  proper  cost  for  the  defense,  and  it 
indeed  uses  similar  code.  The  individual  cell  computation  is  done  in  BDCELL. 
Figure  15  indicates  the  flow  chart  for  BALDEF2. 

This  program  computes  the  balanced  defense  for  a  number  of  costs,  and 
evaluates  each  defense. 

The  defense  computation  is  done  for  a  range  of  lambdas,  for  here  the 
attacker's  lambda  is  the  driving  parameter  for  cost  closing  (see  Appendix  A). 


Subroutine  BDCELL 

BDCELL  is  used  in  the  balanced  defense  procedure  and  finds  the  hardness  and 
cost  for  a  given  cell  according  to  the  balanced  defense  formulation*: 


AREALETH  = 


LAMBDA*  AREA 
'BASE POP 


which  assures  constant  return  for  the  attacker  over  a  number  of  cells.  The  hardness  is 


See  Appendix  C  for  the  formulation. 
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Figure  15.  Flow  Chart  of  BALDEF2 
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BDCELL  (continued) 


computed  from 


and  the  cost  from 


BDPSI  = 


8.1416 


t 


[0.92491  *  (AREALETH50-15-  O.5]3 


BALCOST  =  112.0  +  13.5*  ^ BDPSI 


Saturation  is  checked  at  high  psi  (against  WMIN)  and  low  psi  (against  P5IMIN)  and  the 
highest  psi  is  assigned  in  the  one  case  and  the  unsheltered  posture  in  the  other  (t fie 
unsheltered  posture  is  signaled  with  BDPSI  =  -10),  and  costs  adjusted  accordingly. 
Subroutine  BNDSET 

This  is  the  portion  of  the  bounding  procedure  which  applies  the  upper  and  lower 
bound  theorems.  After  the  program  has  been  through  DEFOPTG  to  find  the  kill,  weapons 
and  cost  for  all  the  lambdas  and  mus,  the  following  are  performed  in  BNDSET: 

1.  For  every  lambda  the  solution  with  cost  closest  to  TCOST  is  selected.  The 
intercepts  are  then  computed,  and  a  printout  is  made  of  this  information. 

2.  To  find  the  upper  bound  the  bracketing  mus  are  used  to  interpolate  an  upper 
bound  at  cost  TCOST. 


3.  Between  each  adjacent  pair  of  lambda  solutions,  it  is  decided  which  lower 
bound  theorem  gives  the  best  bound  —  Pugh's  primary  Pound,  his  secondary 
hound,  or  tire  lower  bound  found  in  this  report.  This  problem  is  merely  a 
geometric  one  and  has  been  reduced  to  tests  on  the  intercepts.  The  lower 
bound  line  is  then  stored. 

4.  From  all  bounding  lines  the  current  lower  envelope  is  computed. 

Subroutine  COLLAPSE 


This  subroutine  collapses  the  hardnesses  available  into  a  new  lethal  area  table 
which  has  no  unused  hardnesses. 


*  From  H.  0.  Horn,  Nuclear  Weapons  Effects:  Approximations  of  the  Phenomena, 
t  An  approximate  analytical  fit  to  the  FINAL  shelters  data. 
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Subroutine  COMFKIL 


COMFKIL  computes 

FKIL  (i,K,J>  =  1.0  -  exp  {  -  W(J,I)  *  D(K)  } 
for  every  hardness  I,  attack  density  K,  and  height  of  burst  J. 

Subroutine  CQRELATE 

This  subroutine  prepares  and  prints  the  following  tables  by  sorting  through  ail 

cells: 

1.  Distribution  of  shelter  mixes  over  the  population  density  brackets. 

2.  For  each  attack  level ,  the  distribution  of  attack  density  over  the  population 
density  brackets.  This  is  done  using  the  lambdas  remaining  from  an  evalua 
tion. 

3.  For  each  attack  level ,  the  distribution  of  attack  density  over  the  various 
shelter  mixes,  again  using  the  lambdas. 

For  samples  of  these  outputs,  see  the  discussion  of  Run  Type  OUTPUT. 

Subroutine  CROSSCOR 

See  Run  Type  CROSSCOR  . 

Subroutine  DCONSTR 

DCONSTR  checks  to  see  if  the  current  defense  for  the  current  cell  is  excluded 
and  sets  DEFEXCL  accordingly.  Normally  it  returns  DEFEXCL  =  0  (not  excluded). 
During  time-phased  runs  it  sets 

DEFEXCL  =  EXCLUDEF  (IDEF,  KDEF) 
where  IDEF  and  KDEF  are  the  current  defense  and  defense  used  for  time-phasing; 
EXCLUDEF  is  computed  in  TIMEFAZE.  If  KDEF  signals  a  split  cell,  DCONSTR 
computes  the  number  of  spaces  of  each  shelter  type  and  checks  defense  option  IDEF 
against  it;  again  returning  DEFEXCL  =  1  (excluded)  or  0  (not  excluded). 

DCONSTR  also  has  entry  ABORT  for  all  improper  terminations  of  the  program. 
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Subroutine-  DEFCOST 


DEFCOST  comixites  the  defense  cost: 

NSTYP 

OCOST  =  £  CSTSHEL  (I)  *  SPACES  (I). 

1  =  i 

Subroutine  DEF1NP 

DEFINE  is  nu  input  routine  which  reads  shelter  parameters,  defense  options , 
and  industry  hardnesses .  The  routine  directly  reads  FDEF,  SPREAD,  SPRIND,  and 
CSTSHEL,  and  computes  COSTPER  for  each  defense  option.  Since  the  input  is  in  a 
difficult  form,  this  routine  is  used  only  for  die  basic  data  decks  and  for  reading  shelter 
data  on  the  punch  decks. 

Subroutine  DEFQPT 

Subroutine  DEFOPT  optimizes  the  defense  for  a  single  cell  for  a  range 
of  current  mu's  (FMUX’s)  and  fills  CDLAG.  the  cell  defense  Lagrangian,  IDFOPT, 
the  optimum  defense  option,  and  CDEFCL,  the  cost  of  the  current  optimum  defense 
option  for  each  FMUX.  It  may  operate  in  one  of  three  modes  -  mixec  lambda,  mini- 
max  lambda,  and  verification. 

In  die  mixed  lambda  mode,  the  task  of  DEFOPT  is  to  choose  one  of  the  15 
defense  options  for  every  FMUX  so  as  to  minimize  the  defense  Lagrangian; 

NAOB  NtAMX  f 

DlAG  WDAOB(JAOB)  WDLAM(I)  *  [PKILLfl)  - 

IAOB-1  1=1 

ALPHA  *  OFFCOST]  +  FMUX(J)  *  DCOST 
where  WOAOB  and  WDLAM  are  defense  weights  on  attack  objective  and  lambda, 
PKILL  and  OFFCOST  are  die  population  kill  and  offense  cost  returned  by  OFFOPT 


112 


DEFQPT  (continued) 

from  the  optimum  attack  for  attack  objective  IAOB  and  lambda  i,  and  DCOST  is  the  cost 
of  the  defense  option  under  consideration.  OFFCOST  is  given  by 

OFFCOST  =  FLAMX(I)  *  D(KDENS)  *  AREA 

where  FLAMX  is  the  lambda,  O(KOENS)  is  the  optimum  attack  density  determined  by 
OFFOPT,  and  AREA  is  the  area  of  the  cell.  ALPHA  is  always  1.0  for  the  double 
Lagrangian  procedure. 

For  each  defense  the  population  kill  and  offense  cost  is  computed  for  an  optimum 
attack.  This  is  done  by  first  calling  GENSPACE,  to  compute  the  number  of  shelter 
spaces  defined  by  the  defense  option;  next  DEFCOST  is  called,  which  calculates  the 
cost  of  the  defense  option  for  that  cell .  The  shelters  are  filled  with  the  population  of 
the  cell  according  to  the  desired  filling  mode  in  FILSPACE,  and  PAYCOEF  then  puts 
the  payoff  (population  and  industry)  into  a  form  easily  usable  by  OFFOPT.  The  optimum 
defense  is  generated  by  OFFOPT  for  a  number  of  lambdas,  and  this  repeated  for  all 
attack  objectives  so  that  the  defense  Lagrangian  in  the  above  form  may  be  computed. 

For  each  mu  DEFOPT  records  the  defense  Lagrangian  CDLAG  and  the  defense 
cost  CDEFCL  of  the  cell  for  the  optimum  defense  chosen  for  that  mu.  Figure  16  shows 
the  fiow  chart  for  DEFOPT  in  the  mixed  lambda  mode. 

In  the  minimax  mode  DEFOPT  operates  as  shown  in  Figure  17.  In  this  mode 
the  object  is  to  consider  several  lambdas  for  the  attacker,  and  to  choose  a  defense  which 
minimizes  the  worst  excess  of  what  is  achieved  for  a  lambda  compared  to  the  best  that 
can  be  achieved.  DEFOPT  first  finds  the  defense  with  the  lowest  Lagrangian  for  every 
lambda  —  it  then  computes  the  delta  or  difference  between  the  Lagrangian  which  each 
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I 


Figure  16.  Flow  Chart  of  DEFOPT,  Mixed  Lambda  Mode 
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DEFOPT  (continued) 


defense  achieves  for  every  lambda  and  the  lowest  Lagrangian  for  that  lambda  --  it  then  finds 
the  highest  delta  for  every  defense  --  and  picks  the  defense  which  minimizes  the  highest 
delta.  The  only  other  computation  is  to  precompute  the  offense  Lagrangian  for  every  lambda 
and  for  every  defense. 

DEFOPT,  in  this  inode,  picks  a  defense  to  minimize 

DELTAMAX  (IDEF)  -  Maximum  {  DELTA  (ILAMX,  IDEF)} 

where 

DELTAflLAMX,  IDEF)  -  DEFLAG  (IDEF)  -  j^i"j.mum  {DEFLAG  (IDEF)} 
and 

DEFLAG  (IDEF)  =  OFFLAG  (ILAMX,  IDEF)  +  FMUX(J)  *  DDCOST  (IDEF) 

NAOB 

OFFLAG  (ILAMX,  IDEF)  =  ^  PKILL  (ILAMX)  -  ALPHA  *  OFFCOST 

IA0B=1 


This  procedure  is  repeated  for  every  mu  until  optimal  defenses  are  found  for  the  cell  for 
every  mu. 

The  verification  mode  of  DEFOPT  is  shown  in  Figure  18.  This  differs  from  the  mixed 
lambda  case  in  that  a  number  of  single  lambda  Lagrangians  are  computed,  and  the  lowest  one 
selected.  The  best  defenses  are  not  saved,  but  only  the  Lagrangian,  cost,  and  weapon 
expenditure.  The  Lagrangian  which  is  minimized  for  every  lambda  and  mu  is 


CURLAG(I,J)  =  Minimum 


|  HLAG(I)  +FMUX(J)*  DCOST} 


where  HLAG  is  computed  by  OF-FQPT  and  is  the  maximum  offense  Lagrangian. 
There  are  no  abnormal  returns  from  DEFOPT. 
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Figure  18.  Flow  Chart  of  DEFOPT,  Verificaticd  Mode 
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Subroutine  DEFOPTG 


DEFOPTG  has  as  its  purpose  to  optimize  the  defense  for  a  group  of 
cells,  indices  ICBEG  and  ICEND,  to  a  total  defense  expenditure  TCOST.  If  it 
is  unable  to  spend  this  amount  (saturation),  it  spends  most  possible  and  exits  with 
ISAT-2  (which  is  otherwise  1)  and  ACOST  equal  to  the  actual  expenditure. 

Actually,  this  is  a  much  too  lofty  statement  of  the  function  performed  by 
DEFOPTG.  In  order  to  optimize  the  defense  for  a  group  of  cells,  DEFOPTG 
requires  the  services  of  an  extremely  able  assistant,  subroutine  DEFOPT.  Recall 
that  on  a  single  cell  DEFOPT  finds  the  optimum  defense  for  a  number  of  mu's 
(or  externally  set  values  for  the  marginal  number  of  people  saved  per  dollar).  It 
is  the  task  of  DEFOPTG  to  drive  DEFOPT  with  a  fixed  set  of  mu's  once  for  each 
ceil  in  the  group  of  cells,  to  accumulate  the  total  cost  of  defenses  optimized  for 
each  mu,  and  to  pick  a  new  set  of  mu's  to  bring  the  achieved  defense  cost  closer 
to  the  target  cost  TCOST.  DEFOPTG  supplies  the  mu's  to  DEFOPT,  whirh  then 
returns  the  cost  of  the  optimal  defense  found  for  each  mu  (other  information  is 
available,  but  DEFOPTG  uses  only  the  cost),  and  DEFOPTG  uses  this  information 
to  pick  a  new  set  of  mu's  after  each  pass  of  optimizing  all  cells.  Thus  we  see  that 
DEFOPTG  contains  only  the  cost  closing  mechanisms  for  adjusting  the  defense  for 
the  group  of  cells  to  a  precise  cost.  Since  individual  cell  optimization  is  done  in 
DEFOPT,  the  various  options  and  strateqems  such  as  expected  attack  levels,  attacker 
objectives,  and  defense  constrains  are  reflected  in  that  subroutine.  Actually 
DEFOPTG  is  a  general  routine,  suitable  for  any  separable  resource  allocation  problem, 
and  in  no  way  considers  the  defense  aspects  of  the  problem. 
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DEFOPTG  (Continued) 

The  procedure  for  cost  closing,  shown  in  the  flow  chart  of  DEFOPTG 
in  Figure  19,  is  to  find  the  optimal  defense  for  each  cell  for  a  number  of  mu 
values.  The  initial  mu's  are  pre-chosen  program  constants  or  input  quantities; 
a  first  pass  is  performed,  optimizing  each  cell  for  this  initial  set  of  mu's.  If 
a  cell  has  the  same  optimum  defense  for  the  entire  mu  range  (this  is  unlikely 
on  the  first  pass,  but  this  is  an  extremely  important  mechanism  on  subsequent 
passes),  its  defense  is  frozen,  and  the  cost  of  its  defense  is  added  to  the  frozen 
cost.  After  the  first  pass  the  two  mu  values  whose  costs  bracket  TCOST  are 
selected  and  a  new  range  of  mu's  is  set  between  those  two  values. 

All  subsequent  passes  follow  the  same  procedure  as  the  first  pass:  each 
unfrozen  cell  is  optimized  for  the  range  of  mu's  and  the  costs  are  noted  for  each  mu. 
If  the  optimized  defense  for  a  cel!  is  the  same  for  the  entire  mu  range,  the  defense 
for  that  cell  is  frozen  and  the  cost  of  the  defense  is  added  to  the  frozen  cost. 

Then  after  the  pass  is  completed  the  bracketing  mu's  are  found  and  the  mu's 
reset  between  them.  On  the  second  and  subsequent  passes,  an  error  quantity  is 
checked  to  see  if  the  procedure  can  be  ended. 

The  desired  outcome  of  the  procedure  is  to  have  an  optimal  defense  which 
is  within  some  small  amount  of  costing  precisely,  TCOST.  However ,  it  may  not 
be  possible  to  get  sufficiently  close  in  by  this  procedure.  It  may  happen  that  for 
very  small  change  in  mu  a  group  of  cells  or  one  sufficiently  large  cell  may  change 
in  cost  enough  to  make  exact  cost  closing  possible.  Thus  it  is  reasonable  that  the 


error  criterion  is 

ERRZ 


PLTMU  *  DELTC  4  Q  0Q01 
FMULOW  TCOST 
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CALL  DEFOPTG 

Group  defense  optimization  1 

procedure  (DEFOPTG) - 1 - 

i  Set  initial  m 's  Clear  Defense 
I  Strat.  Clear  iteration  ctr.  Sel 


t»get  cost  TCOST.  Find  lowest 
defense  A  for  use  in  OFFOPT 


Figure  19.  Flow  Chart  of  DEFOPTG 
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Figure  19  (continued).  Flow  Chart  of  DEFOPTG 
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DEFOPTG  (Continued) 

where  DLTMU  is  the  change  in  mu  between  upper  and  lower  values,  and 
FMULOW  is  the  lower  mu  value,  DELTC  is  the  fractional  change  in  mu  and 
DELTC/TCOST  is  the  fractional  change  in  cost  for  the  two  values.  Thus  if 
either  fractional  quantity  is  small,  i.e.,  if  either  mu  closing  or  cost  closing  is 
close,  then  the  iterating  will  cease.* 

When  the  error  criterion  has  been  satisfied  (but  the  cost  is  still  not 
closed  exactly),  two  mu's  have  been  found,  one  with  too  high  a  cost  and  one 
with  too  low  a  cost.  Most  of  the  cells  will  have  the  same  defense  for  both  mu's, 
and  will  have  been  forzen.  Since  the  cells  which  have  not  been  frozen  yet  are 
flipped,  one  by  one,  from  the  lower  cost  defense  to  the  higher  cost  defense,  and 
the  additional  cost  for  that  cell  is  added  to  the  amount  spend  in  the  lower  cost 
defense.  This  process  of  flipping  cells  to  the  higher  defense  is  continued  until 
the  cell  is  reached  where  expended  cost  exceeds  TCOST.  Then  this  cell  is  split, 
putting  part  at  the  higher  defense  andpart  at  the  lower  defense.  Finally,  all 
remaining  cells  are  set  to  the  lower  cost  defense  (again,  between  choices  generated 
by  the  two  bracketing  mu's),  and  the  defense  is  set  for  the  entire  group  of  cells. 

If  after  the  first  pass,  the  costs  computed  for  the  initial  mu  range  do  not 
encompass  TCOST,  one  of  two  things  happen--if  the  computed  costs  are  too  low, 
the  saturation  indicator  (I SAT)  is  set  and  all  cells  are  set  to  the  defense  with  the 
highest  computed  cost.  If  the  computed  costs  are  too  high,  the  mu  range  is  expanded 

*ERRZ  has  a  more  theoretical  justification:  DLTMU*  DELTC  is  an  upper 
bound  on  possible  difference  is  fatalities  between  the  two  points  generated  by  the 
two  mu's  in  question  (from  the  definition  of  mu),  and  FMULOW*  TCOST  is  a  lower 
bound  on  the  number  of  people  saved  ay  a  shelter  program  of  cost  TCOST.  Hence 
ERRZ  is  an  upper  bound  on  the  fractional  error  in  fatalities  from  using  either  one  of 
the  mu's. 
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DEFOPTG  (Continued) 


by  10  and  the  first  pass  is  repeated. 

For  the  operation  of  DEFOPTG  in  the  verification  procedure,  see  Figure  20  . 

In  this  procedure  the  idea  is  to  find  optimum  defenses  for  a  range  of  20  lambda's 
and  20  mu's,  and  to  save  the  total  kill,  weapon  expenditure,  and  cost  expenditure 
for  each  of  the  400  lambda-mu  combinations.  For  verification,  exact  cost 
closing  is  not  required;  hence  this  mode  of  operation  is  not  iterative  in  DEFOPTG; 
the  subroutine  makes  one  pass  through  all  the  cells,  accumulating  totals  for  each 
lambda-mu  combination,  and  that  is  all  that  is  required  for  verification. 

There  are  only  two  abnormal  returns  from  DEFOPTG.  One  is  in  the 
verification  procedure--if  more  th'  .  20  lambdas  or  mus  are  used,  the  program  aborts. 
The  other  is  shown  at  the  very  end  of  Figure  20  :  if,  because  of  rounding  perhaps, 
ail  the  cells  are  flipped  and  set  to  the  upper  defense  without  finding  the  one  to  split, 
then  if  the  upper  defense  cost  is  close  to  TCOST  then  the  routine  returns  normally 
(with  all  cells  flipped  to  the  higher  defense).  If  the  upper  defense  cost  is  not  close, 
the  routine  aborts. 

Subroutine  DEFPRNT 

This  is  a  print  routine  responsible  for  printing  the  defense  when  it  is  input;  for 
a  sample  print  see  the  discussion  of  Run  Type  INPUT. 

Subroutine  DEFPUNCH 

This  subroutine  punches  a  defense  onto  cards  for  saving.  First,  a  header  is 
punched,  and  then  the  defense  options,  40  to  a  card.  Next,  split  cell  information,  and 
all  parameters  and  shelter  data  are  punched.  For  a  list  of  the  parameters  punched  see 
DEFSPEC  in  the  data  sequence  definition  section. 
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Optimize  defense 
Vm  (in  DEFOPT) 
Accumulate  costsVu 


Figure  20.  DEFOPTG  in  Verification  Mode 


Subroutine  DEFSPEC 


See  Run  Type  DEFSPEC  for  functions  of  this  subroutine. 

Subroutine  EVAL 

The  task  of  EVAL  is  to  find  the  optimum  attack  against  all  cells  for  a  range  of 
attack  levels  and  compute  the  attack  payoff,  population  kill,  and  industry  kill  for  each 
attack  level.  Although  EVAL  is  used  tc  evaluate  a  defense  for  different  times  of  attack 
(IAPQP),  attack  objectives  (!A0B),  and  filling  modes  (IFILL),  these  parameters  are  not 
reflected  in  EVAL  but  rather  in  FILSPACE  and  PAYCOEF,  which  are  called  from 
OFFOPTZ .  The  flow  chart  of  EVAL  is  shown  in  Figure  21 . 

The  portion  of  the  evaluation  task  performed  in  EVAL  is  to  select  lambdas  to 
drive  OFFOPTZ,  finally  selecting  lambdas  which  give  weapon  expenditures  close  enough 
to  the  desired  attack  levels.  For  the  initial  preset  set  of  lambdas  OFFOPTZ  finds  the 
optimum  offense,  returning  weapon  expenditures  for  each  lambda.  If  the  range  of  expendi¬ 
tures  does  not  encompass  the  desired  attack  level  range,  the  lambda  range  is  increased 
and  the  initial  pass  is  redone.  When  the  weapon  expenditure  range  encompasses  the 
attack  level  range,  the  procedure  is  simply  to  find  the  lambdas  with  expenditures  bracketing 
each  desired  attack  level,  setting  lambdas  linearly  between  the  bracketing  lambdas,  and 
finding  a  new  set  of  optimum  attacks  in  OFFOPTZ.  After  each  pass  the  kill  at  each 
precise  attack  level  is  estimated  from  the  bracketing  values  for  each,  and  an  error  bound 
is  determined.  If  the  error  is  small  enough,  the  procedure  terminates;  otherwise,  new 
lambdas  are  set  linearly  between  the  bracketing  lambdas  and  a  new  set  of  optimum  attacks 
is  found  in  OFFOPTZ,  until  the  error  is  small  enough. 
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Evaluation  of  defense  for  range  of  CALL  EVAL 

attack  levels  (EVAL)  c - - 

i — * — - - ™  ■  -sfeet  Initial  x's 


Figure  21.  Flow  Chart  of  EVAL 


EVAL  (continued) 


Ai  Ae 

Figure  22.  Estimation  and  Error  for  Kill  at  Desired  Attack  Level 

Figure  22  shows  the  estimation  procedure  and  error  criterion  for  EVAL. 

F^  (A^)  and  ^2^2^  are  tvv0  level s  corresponding  to  lambdas  which  bracketed 

the  desired  attack  level  Ae.  The  estimate  used  for  fatalities  at  Ae  is  the  linear  inter¬ 
polation  between  F^  and  F^.  The  point  Fe  is  a  lower  bound  on  fatalities  at  Ae  since 
the  payoff  curve  is  convex  .  The  lines  upper  bound  #1  and  upper  bound  #2  have  slopes 
nf  the  bracketing  lambdas  (see  the  section  on  double  Lagrange  theory  for  why  they  are 
upper  bounds),  and  these  lines  provide  an  upper  bound  for  the  payoff  at  any  attack  level. 
To  provide  an  upper  bound  for  the  payoff  at  Ae,  we  simply  take  the  lower  of  the  two 
upper  bounds.  The  possible  error  is  thus  as  indicated,  and  this  error  is  converted  into 
an  error  in  attack  level  using  the  lambda  value  of  the  upper  bound. 

X”— — — —  —  - - — — "  - — — 

Convex  because  of  decreasing  returns  and  because  the  attacker  could  split  attack 
among  cells  to  achieve  any  intermediate  point. 
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EVAL  (Continued) 


The  decision  to  terminate  the  iteration  is  made  on  that  quantity.  For  information  purposes, 
the  lambda  value  at  Ae  is  estimated  as  (F2  -  F  i>/(A  2  "A]). 

There  is  only  one  abnormal  return  from  EVAL  and  it  can  occur,  as  it  turns  out, 
if  on  the  first  pass  the  lowest  attack  level  is  precisely  equal  (to  the  last  bit)  to  the 
iowest  weapon  expenditure. 

Subroutine  FILS  PACE 

FILSPACE  first  determines  the  number  of  people  in  the  cell  at  time  of  attack 
(from  IAPOP)  and  then  fills  the  shelter  spaces  in  one  of  three  modes  —  IFILL  =  1, 

2,  or  4  —  optimal,  fractional  alert,  and  balanced  defense.  The  routine  begins  with 
the  SPACES  array  and  fills  the  population  into  the  POPSH  array. 

Optimal  filling  mode  fills  the  population  into  the  hardest  shelters  first,  working 
down  until  either  population  or  space  is  exhausted.  Any  remaining  population  is  placed 
in  the  unsheltered  category. 

The  fractional  alert  mode  assumes  that  some  fraction  of  the  population  has  received 
the  alarm.  If  the  alerted  population  is  more  than  available  space,  all  spaces  are  filled. 

If  the  alerted  population  is  less  than  available  space,  each  shelter  is  filled  the  same 
fraction  of  capabity.  All  unalerted  population,  of  course,  is  unsheltered. 

The  balanced  defense  mode  is  not  really  a  filling  mode,  but  rather  a  way  of 
implementing  the  variable  hardness  requirement  for  balanced  defense  in  the  evaluation 
process .  If  the  cell  is  unsheltered ,  operation  of  FILSPACE  is  normal .  If  not,  subroutine 
BOCELL  is  called  which  computes  the  proper  hardness  and  its  lethal  area.  FILSPACE 
then  sets  the  proper  kill  fractions  by  calling COMFKIL ,  and  operation  is  then  normal. 
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Subroutine  GENSPACE 


GENSPACE  computes  the  number  of  shelter  space*  in  a  ceil  indicated  by 
each  defense  option.  It  first  finds  the  base  population  BASEPOP  using  IDPOP  and 
then  computes  each  shelter  space  from 

SPACES(I)  =  FDEFU,  IDEF)  *  BASEPOP. 

Subroutine  GETCELL 

This  subroutine  reads  ceil  data  in  the  form  shown  in  option 
LIVE  DAT  A. 

Subroutine  LAMSET 

LAMSET  uses  FLAM(I),  and  FLAM(NLAM)  and  computes  intermediate 
FLAM's  in  an  ascending  geometric  series.  LAMSET  also  contains  a  diagnostic  print 
routine  as  entry  LOBNDZ , 

Subroutine  MIXLAM 

See  Run  Type  MIXLAM  for  this  subroutine. 

Subroutine  Of  --,l0T 

OFTOPT  generates  the  optimum  attack  density,  height  of  burst, 
population  kill,  and  industry  kill  on  a  single  cell  for  a  set  of  NLAMX  attack 
Lagrange  multipliers  (FLAMX's).  The  flow  chart  is  shown  in  Figure  23. 

The  basic  task  of  OFFOPT  is  to  find  the  attack  density  index  and  the 
height  of  burst  index  which  maximizes  the  offense  Lagrangian 
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Figure  23.  Flow  Chart  of  OFFOPT 


OFFOPT  Com i rued) 


NCOMP 

£  A(L)  *  FKIL  [iCOMP  (L),  KDENS,  JHOb] 

L  -  1 

-  FLAMX  (ILAM)  *  AREA  *  D(KDENS) 

for  every  Lagrange  multiplier  FLAMX  (ILAM).  See  the  glossary  for  definition 
of  terms.  For  every  optimum  attack  so  found,  OFFOPT  saves  the  optimal 
attack  density  IDOPT,  height  of  burst  IHQBOPT,  payoff  HQPT,  offense 
Lagrangian  HLAG,  population  kill  PKILL,  industry  kil*  VKILL,  and  weapon 
expenditure  WEPEQ. 

Three  time  saving  short  cuts  were  added  to  this  procedure:  1)  FKIL 
is  not  used  as  a  three-dimensional  array  but  rather  as  15  two-dimensional 
arrays;  2)  in  the  attack  density  scan,  if  the  optimum  height  of  burst  is  ground 
burst,  all  subsequent  higher  densities  are  ground  burst  also;  3)  in  the  attack 
density  scan,  if  for  the  lowest  (and  hence  all)  lambda  the  Lagrangian  reverses 
and  starts  decreasing,  all  subsequent  higher  densities  will  have  lower  Lagrangians, 
and  the  scan  is  terminated. 

There  are  no  abnormal  returns  from  OFFOPT. 

Subroutine  OFFOPTZ 

OFFOPTZ  finds  the  optimum  attack  for  all  cells  for  a  range  of  externally  set 

lambdas  (FLAMX).  OFFOPTZ  is  a  simple  routine  and  fora  normal  (unsplit)  cell,  the 

procedure  is  to  read  the  cell,  and  next 

CALL  GENSPACE 

CALL  FILSPACE 

CALL  PAYCOEF 

CALL  OFFOPT 


Maximum 

JHOB 
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OFFOPTZ  (Continue!) 

and  lastly  to  accumulate  the  total  attacker's  payoff,  population  kill,  industry  kill, 
and  weapons  expended  for  every  lambda.  All  the  (layoffs  and  kills  are  returned  from 
subroutine  OFFOPT 

All  cells  are  processed  in  the  above  manner  except  for  the  split  cells,  of  which 
there  is  ususally  one  in  a  defense.  (The  split  cel!  is  generated  by  DEFOPlG  for  final 
cost  closing.)  OFFOPTZ  processes  that  cell  twice,  weighting  the  kills  for  the  primary 
and  secondary  defense  by  the  fraction  of  the  cell  at  each  defense.  Figure  24  shows  the 
flow  chart  for  the  entire  procedure. 

There  are  two  abnormal  returns  from  OFFOPTZ:  if  a  cell  has  defense  option 
zero,  the  program  aborts;  and  if  a  cell  has  a  negative  defense  option  (indicating  a  split 
cell),  but  the  cell  is  not  contained  in  the  split  celi  list,  the  program  aborts. 

Subroutine  PAYCOEF 

PAYCOEF  first  computes  the  industrial  value  and  population  in  each  hardness 
component  (using  POPSH  computed  in  FILSPACE)  and  stores  these  in  array  VALINO 
and  POPUL.  It  then  computes  the  value  at  each  hardness  component  for  the  current 
attack  objective  using  POPMULT  and  VINOMULT,  and  stores  the  attacker's  value 
at  each  hardness  in  a  compacted  array  A  (having  no  zero  components),  using  an  array 
ICOWIP  to  store  the  hardness  components  represented  in  A. 
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CALL  OFFOPTZ 


Figure  24 .  Flow  Chart  of  OFFOPTZ 
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Subroutine  PERTURB 


The  purpose  of  this  subroutine  is  to  adjust  a  i  lethal  areas  in  response  to  a 
fractional  change  in  hardness  of  DELTA.  The  equation  used  is 


W(J,i)  =  WU,I)  -  DELTA*  PSKI)  * 


2  *  W(J  ,\) 

BETA  *  PSKI)  -  X 


I  x  -  .00001875  *  HQB(J)  *  ALPHA 

W*re  [0.34303*  yjmj]  BETA-0.75 

where  ALPHA  is  8.7  or  10.0  and  BETA  is  2.28  or  2.0,  depending  on  whether  the 
lethal  area  is  greater  or  less  than  3.03.  The  above  equation  is  found  by  taking  partial 
derivatives  from 


PS |  =  ALfii.A-  [l.0+  H0B2  +  A(1  -  .25  D0-75)]  * 

□BETA  L  J 

where  here  A  and  D  are  height  of  burst  and  distance  in  kilofeet  for  a  1  kt  weapon. 

Subroutine  PICKLAM 

This  subroutine  is  responsible  for  choosing  the  lambdas  for  the  bounding  procedure. 
The  procedure  used  is  exactly  that  described  in  the  theoretical  section  on  the  bounding 
procedure.  The  lambda  sweep  begins  in  the  middle  and  extends  both  directions;  hence 
two  sets  of  10  lambdas  are  chosen  sweeping  In  opposite  directions  from  the  previous  ends 
of  the  sweep. 

The  controlling  parameter  for  the  separation  between  the  lambdas  is  EPS;  it 
enters  strongly  into  the  expressions  for  the  new  lambdas. 


Ifrom  H.  O.  Norn,  'I'JucreaFWeaponsEffects:  Approximations  of  the  Phenomena" . 
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Suijr  oulme  POPQUANT 

This  subroutine  generates  nine  population  density  brackets  and  sorts  the  cells 
into  the  nine  brackets ,  accumulating  totals,  it  is  a  prelude  to  subroutine  CORELATE; 
an  example  of  the  printout  is  the  first  block  of  print  in  Run  Type  OUTPUT. 

Subroutine  PREBND 

This  is  the  driver  for  the  bounding  procedure.  From  the  flow  chart  in  Figure  25, 
it  can  be  seen  that  the  basic  cycle  is 

1.  Select  lambdas  in  PICKLAM 

2.  Find  Lagrangian  solutions  in  DEFOPTG 

3.  Compute  bounds  in  BNDSET 

4.  Repeat 

The  variations  on  this  theme  are  in  the  different  modes;  these  modes  are  described  in 
the  section  on  run  type  BOUND . 

Subroutine  PRNTBLK 

PRNTBLK  evokes  the  entire  range  of  diagnostic  block  print  routines  from 
ATTOBJZ  to  SPLITDAZ. 

Subroutine  PRTOCT 

PRTOCT  is  an  octal  print  routine  with  entries  PR  TINT  for  integer  print 


and  PRTDEC  for  decimal  print. 


Figure  25.  Flow  CKart  of  PREBND 
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Function  RANDNOO 


RANDN  is  a  pseudo-random  number  generator. 

Subroutine  RDATT 

RDATT  is  the  input  routine  which  reads  heights  of  burst,  hardnesses  and 
either  lethal  areas  or  lethal  radii.  It  then  reads  the  actual  heights  of  burst  to  be 
scanned  and  prints  both  the  entire  lethal  area  table  and  a  reduced  table  with  heights 
of  burst  to  be  used.  The  input  format  is  included  in  the  INPUT  data  sequence  definition. 

Subroutine  RDCELDT 

RDCELDT  is  used  in  reading  the  basic  data  deck,  and  it  reads  variables 
FINDMAX ,  AREAMAX ,  POPMAX ,  and  then  reads  the  cell  data  for  each  cell . 

Subroutine  RDCELL 

RDCELL  unpacks  the  basic  ell  data  for  cell  ICELL  from  storage  in  NCELA 
and  NCELB  and  stores  it  in  usable  form  in  common  block  CURCELL. 

Subroutine  RDCTDAT 

This  subroutine  reads  trac'  data  in  the  form  shown  in  option 
LVCTDAT. 

Subroutine  SHELINP 

This  subroutine  reads  in  shelter  data  in  the  form  specified  in  the  data  sequence 
for  R^n  Type  INPUT. 

Subroutine  5TATCQST 

See  Run  Type  STATCOST  for  a  discussion. 
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Subroutine  SPLITDAZ 


SPLITDAZ  contains  20  of  the  diagnostic  block  print  routines,  each  under  a 
separate  entry  point. 

Subroutine  STCELL 

STCELL  packs  the  basic  cell  data  for  cell  ICELL  into  NCELA  and  NCELB 
from  CURCELL  in  a  form  usable  by  RDCELL. 

Subroutine  STCELLDT 

This  subroutine  tests  the  population  and  industrial  density  cf  cells 
read  in  GETCELL  and  either  stores  the  ceil  or  accumulates  it  into  tail 
cells. 

Subroutine  TABLDENS 

The  function  of  TABLDENS  is  described  in  that  Run  Type  discus,.... .  The 
p/a  breakpoint  between  two  attack  densities  is 

p/x  =  PQ)  -  DC2) _ 

ATTPAY1  -  ATTPAY2 

Where  the  D's  are  the  attack  densities  and  the  ATTPAY's  are  the  fractional  kills. 

Subroutine  TABLE 

This  is  a  table  lookup  subroutine. 

Subroutine  TIMEFAZE 

Discussed  under  Run  Type  TIMEFAZE. 

Subroutine  WORKOUT 

For  discussion  see  Run  Type  EXERCISE. 


138 


1 


UNCLASSIFIED 


Vrunty  Cliwftcitiftw 


DOCUMENT  CONTROL  DATA  RAD 


LAMBDA  Corporation 


AN  OPTIMIZATION  STUDY  OF  BLAST  SHELTER  DEPLOYMENT,  VOLUME  III: 
APPENDIX  H  =  "BLAST"  -  THE  COMPUTER  PROGRAM 


I*  Otic  Rip  Ti  »L  mo  Tfcl<  l>p*  at  —pan  mad  )ne  hr*  rrr  t fmtrwj 


la  *u  tmoRi3i  tFitn  n#m«,  **<<3<fl+  inUlm J,  iamt  m«<«) 


David  L.  Mitchell 


1«  RtPOR  T  DA  Tit 


l»«.  total  no  or  PACti  It*,  no  or  rifi 


■  A  CONTRACT  on  &RANT  no 

OCD-PS-66-113 


Work  Unit  1632A 
Subcontract  138-5 


IN.  ORiQ<N«TOR‘l  REPORT  huk«ih|I| 


Report  3 


O  ▼  M r  «  REPORT  ND<»  (Ant  *+tat  iwNr*  Mr  A*r  M  —mlgnad 
thtu  t+port) 


rTr  o.ji^.BuTioH  it.tem.ht  document  ha s  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited. 


M  iUPPHUlNTARr  NOm 


1*1  (PORIDRINC  IflUTIRT  ACTIVITV 


Office  of  Civil  Defense 
Department  of  the  Army 

"'"‘"This  study  examines  methods  of  determining  blast  shelter  de¬ 
ployments  and  of  assessing  their  performance  for  a  variety  of  nuclear 
attacks.  The  goal  is  not  to  seek  a  single  optimal  deployment,  which 
generally  requires  making  arbitrary  assumptions  on  the  nature  and 
size  of  the  attack,  thus  overlooking  the  attacker's  freedom  of  choice 
after  a  blast  shelter  program  has  been  deployed.  Rather,  the  study 
seeks  "stabilized"  deployments  which  protect  population  almost  as 
well  as  an  optimal  deployment,  even  though  it  is  not  truly  optimal 
for  any  specified  attack. 

The  study  examines  the  attacker's  freedom  to  vary  force  level, 
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